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ABSTRACT 

We have used the SINFONI integral field spectrograph to map the near-infrared K- 
band emission lines of molecular and ionised hydrogen in the central regions of two 
cool core galaxy clusters, Abell 2597 and Sersic 159-03. Gas is detected out to 20 kpc 
from the nuclei of the brightest cluster galaxies and found to be distributed in clumps 
and filaments around it. The ionised and molecular gas phases trace each other closely 
in extent and dynamical state. Both gas phases show signs of interaction with the 
active nucleus. Within the nuclear regions the kinetic luminosity of this gas is found 
to be somewhat smaller than the current radio luminosity. Outside the nuclear region 
the gas has a low velocity dispersion and shows smooth velocity gradients. There is no 
strong correlation between the intensity of the molecular and ionised gas emission and 
either the radio or X-ray emission. The molecular gas in Abell 2597 and Sersic 159-03 
is well described by a gas in local thermal equilibrium (LTE) with a single excitation 
temperature Tcxc ^ 2300 K. The emission line ratios do not vary strongly as function 
of position, with the exception of the nuclear regions where the ionised to molecular 
gas ratio is found decrease. These constant line ratios imply a single source of heating 
and excitation for both gas phases. 

Key v^rords: galaxies:clusters:individual: Abell 2597 - galaxies: clusters: individual: 
Sersic 159-03 - cooling flows - infrared: ISM. 



1 INTRODUCTION 

Cool cores are regions at the centre of rich clusters where 
the hot thermal X-ray emitting gas ( T ~ 10* K) is dense 
enough to cool radiatively within a Hub ble time (see 
iPeterson fc Fabianl [200^ : iFabian et all 1 1994 for reviews). 
Cooling rates of the order of 100 M© yr~^ and up to 
1000 M0 yr^'^ have been estimated for this hot X-ray 
gas (e.g., iPeres et al.lll99§ ). However, recent Chandra and 
XMM-Newton X-ray spectra show that little or no X-ray 
emitting gas (<10 %) cools below one third of the virial 
temp erature (e.g., iKaastra et al.l I2OOII : iPeterson fc Fabianl 
I2OO6I ). The solution most often invoked in the literature is 
that some form of reheating balances the radiative cooling 
of the X-ray gas. 



be explained by the simple cooling of the intracluster gas 
through this temperature i.e., additional heating is needed 
ijFabian et al.lll98ll : iHeckman eral]| 19891 ) . 

More recently, work in the infrared and mm- 
wavelengths has shown that there are large quanti- 
ties of molecular gas at the centre of these clusters 
with temperatures between 15 and 2500 K extending 
at l east up to 20 kpc f rom the BCG nuclei (e.g., 
J05; Jaffc fc Bremeij 1 19 97': Jaffc. Breme r fc van der Werj 
2001; Falckc et all Il998 : Edge 2001; Edge et all I2OO2I 



Wilnian et al.l |2002|: ISalome fc Combesi l20of: iHatch et all 



Johnstone et al. I l2007l : IWilman. Edge fc Swinbaii^ 



years 
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Substantial cooler gas and dust com ponen ts 
in the cores of these gal axy clusters (|Edg el 



[ ;as has a cooling time of order 
19831 ; iBlack fc van Dishoeckl [l987l : 



2001 



2002 



Irwin, Stil fc Bridge^ I2OOII : ISalome fc Combesi 
O'Dca et al. 200^). Extended 10* K emission-line nebulae 
are found surro u nding Brightest Cluster Gal a xies (BCG) 
llHeckman et al.l Il989l ; ICrawford et al. 1 Il999l ; iJaffe et al.l 
l2005l . herafter J05). These nebulae are observed to extend 
at least up to 50 kpc from the BCG nuclei (J05). This 



iMalonev" et al.lll996l ). Without some form of heating one 
would expect this gas to collapse rapidly and form stars. 
Although there is strong evidence that starformation does 
take place at the centres of cool core galaxy clusters it 
is not yet observed to match the extended ga s nebulae 
ijMcNamara fc 0'Connelllll992l ; lO'Dea et al.|[200a Oonk et 
al. in prep.). 

The heating and cooling of the molecular and ionised 
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gas phases are important elements in the energetics of the 
cool core region. An energy source comparable to that 
needed to stop the hot X-ray gas from cooling is necessary to 
heat these colder phases (J05). The primary source of ionisa- 
tion and heating of the H2 and H II must be local to the gas 
(JOS; I Johnstone fc FabianI Hoisl l. consistent with a stellar 
photoionising source. However, stars are unable to explain 
the h igh temperature of the ionised gas (|Voit fc Donahug 
I1997I . hereafter VD97). The molecular hydro gen lines are 
much stronger relative to the ionised hydrogen lines than in 
other types of extragala ctic sources, such as AG N or star- 
burst galaxies (e.g., J05: lDavies et al.ll2003l . l2005h. 

The ratio of H2 to HII emission lines (JOS; iHatch et al.l 

as well as detailed analysis of the mid-inf rared and 
optical line ratios (VD97; Ijohnstone et al.l l2007h indicate 
that to explain both heating and ionisation balance, photons 
harder than those available from O-stars are needed. How- 
ever, often very high ionisation lines are missing (e.g., [Ne V] 
- typical of AGN spectra. If present, the source of these pho- 
tons is elusive. Alternatively shock heating has been consid- 
ered, however the characteristic [O HI] 4363 Angstrom line 
is missing (VD97). High energy particles hav e been invoked 
to explain this problem l|Ferland et al.l [20091 ). It is of great 
importance to pinpoint the nature of the heating mechanism 
and include it in models of cooling flows in galaxy clusters 
as well as models of galaxy growth and evolution. 

Cool core clusters are the low redshift cluster scale 
analogs of high redshift galaxy scale cooling flows. To un- 
derstand the formation of massive galaxies at high redshift 
and the feeding and feedback mechanisms in AGNs it is im- 
portant to understand the heating of the cool gas in BCGs. 

All gas phases observed in the intracluster medium 
require reheating to avoid catastrophic cooling. A vari- 
ety of heat sources to counteract this cooling gas have 
been proposed o ver t he years: AGN feedback (e.g.. 



Churazov et all I2OO1I; 



Blanton. Sarazin fc McNamaral 



20031: iDaUa Vecchia et~aLl |2004 iMcNamara et all I2OOII: 



Birzan et al.l 2004 ) , low velocity shock waves ( Fabian et al.l 
20061 ), co nduction (VD97), h ot stars (VD97) and energetic 



particles l|Ferland et al.|[2009l ). None of these heat sources 
have so far been able to match the detailed character- 
istics of cool core galaxy clusters. Whatever the heat 
mechanism, the cooling region extends over hundreds of 
kpc across the cluster core, and heating is unlikely to 
balance the cooling exactly over such a large region. Some 
residual cooling will occur and presumably corresponds 
to the emission line nebulae and star clusters surrounding 
brightest cluster gal axies (BCG) at the centres of cool cor e 
galaxy clusters (J05: lHatch et al.ll2005l : rRaffertv et a"l]|2006l ). 

The BCGs at centres of c ooling clusters fall within a 
regio n of the BPT dia g ram (|Baldwin. PhiUips fc TerlevichI 
I2OO4I : ICrawford et~aLl Il999l . VD97) that is occupied 
by LINERs and AGN. In our previ ous samples (JOS; 
iJaffe k. Bremed Il997l : Ijaffe et all I2001I ) we have focused 
on the LINER-like BCGs and we continue to do so here. 
These clusters were originally selected based on their high 
cooling rates, strong Ha, H2 emission and low ionisation 
radiation. LINER-like BCGs were chosen because we wish 
to minimise the role that their AGN have on the global 
radiation field. In the work presented here we focus on two 
LINER-like BCGs from our previous samples, PGC 071390 



in Abell 2S97 (hereafter A2597) and ESQ 0291-G009 in 
Sersic 1S9 -03 (h e reafte r S1S9). Optical ([O III]/H^ and 
[O I]/Hq, iBakerl (|200Sl . VD97;)) and mid-infrared spectra 
([Ne III]/[Ne II] and [Ne V]/[Ne II], Jaffe & Bremer in 
prep.) of these BCGs indicate that their ionising sp ectrum 
is ve ry soft i.e. they are extreme LINERs (VD97; iBakerl 
l2005h . 

However, these BCGs do contain radio-loud AGN. 
Their 1.4 GHz radio specific luminosity is, 2 9.3x10^^^ and 
1.6x10^^ erg s'^ Hz~^ for A2S97 and S1S9 (|Birzan et all 



I2OO8I ) respective ly, which are typica l for BCGs in cool core 
galaxy clusters ( Quillen et al.l |2008| ) . In this work we will 



concern ourselves with the extended molecular and ionised 
gas surrounding the BCGs in A2S97 and S1S9. 

A2S97 and to a lesser extent S1S9 have been the sub- 
ject of numerous investigations and have been observed at 
many wavelengths frorn radio to X-rays (e .g., JOS; VD97; 
lO'Dea et al.lll994 |2004| : IClarke et al.H200Sl ). In both clus- 
ters ionised and molecular gas was observed to at least 
SO kpc and 20 kpc fro m their BCG nuclei respectively (JOS; 
iHeckman et al.lll989l ). Previous investigations of these ob- 
jects made use of narrowband imaging and longslit spectra. 
Using long slit observations we were only able to sample 
parts of the extended line emission and with narrowband 
imaging no information on the dynamics of the gas is ob- 
tained. Furthermore the emission sampled with longslits in 
previous observations is often strongly dominated by strong 
emission from the BCG nucleus. As we will show below, the 
emission away from the nuclei has very different character- 
istics. 

There are a number of kinematic problems con- 
cerning the cooler gas phases in cool core clusters. 
In nearby clusters ionised and molecular gas is found 
in thin, long lived, multi-kpc scale, filamentary struc- 
tures surrounding the BCG (e.g., Fabian et al.l I2OO8I : 



iHatch et all l2005l : ICrawford erall |2005| ). These structures 
show smooth velocity gradients, but no rotation beyond 
the central few kpc (Wilman. Edge fc Swinban^ l2006l : 
iHatch. Crawford fc FabianI l2007l ). The molecular clouds in 
these structures are dense and without kinematic support 
should fall to the galaxy centre. However, they show no 
signs of infalling. Velocities barely exceed 100 km s~^ (JOS; 



^ ely 

IWilman. Edge fc SwinbankI |2009| . and this work), whereas 
infall velocities should ex ceed ^600 km s"'^ . Magnetic sup- 
port has been invoked bv i Fabian et al.l (|2008| ) but there is no 
observational evidence yet for the required ordered magnetic 
fields in clusters. There has also been speculation whether 
or not all the molecular and ionised gas is locked up in these 
dense filaments or if a more diffuse underlying component 
exists. 



1.1 This project 

Here we present the first deep K-band integral field (IFU) 
spectroscopic observations of the cluster cores in A2S97 and 
S1S9, taken with the Spectrograph for INtegral Field Obser- 
vations in the Near-Infrared (SINFONI) on the Very Large 
Telescope (VLT). With these observations we sample the 
molecular and ionised gas phases over a major fraction of 
each cluster's BCG. Our observations are able to provide in- 
formation on the distribution, kinematics and temperature 
of this gas. Using these measurements we can compare in de- 
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tail the kinematic and thermal structure of the gas with the 
X-ray and radio structures, which represent respectively the 
primary source of mass in the environment and the primary 
source of local energy input. Similar data on t hree other cool 
core c lusters has recently been presented by IWilman et alj 
but we are the first to make a detailed comparison 
of such data with radio and X-ray observations of cool core 
clusters 

In Section 2 we describe the observations and the data 
reduction. We discuss the morphology and kinematics of the 
molecular and ionised gas in A2597 in Sections 3 and 4 and 
similarly for S159 in Sections 5 and 6. In Section 7 we will 
discuss the excitation of the molecular gas and in Section 

8 we compare the observed gas structure to high resolution 
X-ray and Radio maps. We summarise our results in Section 

9 and present our conclusions in Section 10. 

Throughout this paper we will assume the following 
cosmology; Ho—72 km s~^ Mpc~^, r2in=0.3 and r2A=0.7. 
For Abell 2597 at z=0.0821 (VD97) this gives a luminosity 
distance 363 Mpc and angular size scale 1.5 kp c arcsec"^. 
For Sersicl59-03 at z=0.0564 ijMaia et al.lll987l 'l this gives 
a luminosity distance 245 Mpc and angular size scale 
1.1 kpc arcsec^^. 



2 OBSERVATIONS AND REDUCTION 
2.1 Near Infrared Data 

The near infrared (NIR) observations were performed in 
the K-band with the in t egral field spectrogr aph SINFONI 
lEisenhauer et al.l l2003l : iBonnet et al] |2004| '1 on the Very 
Large Telescope (VLT). SINFONI is an image slicer, slic- 
ing the image into strips before dispersing the light using 
32 slitlets. The instrument has a spectral resolving power 
oi R K. 4000 in the K-band. Opting for a 8"x8" field 
of view (FOV) the spatial pixels each cover an area of 
0.125" xO.250". Each spectral pixel covers 2.45x10"* ^.m 
in wavelength, oversampling the resolution by a factor two 
(i.e. Nyquist sampling). The total on-source exposure time 
for each source is listed in Table [T] 

The observations were done in a 'ABBA' pattern 
(A=source, B=sky) and each set was followed by a pointing 
observation to keep track of pointing drifts. All observations 
were done such that the FOV was oriented with north point- 
ing up. The spatial extent of each slitlet is then oriented east- 
west. Equal amounts of time were spend on the sky and on 
the source. Each science observation has an exposure time 
of 600 seconds. Each pointing observation has an exposure 
time of 60 seconds. The observations were carried out in 
July and August of 2005 in photometric sky conditions with 
a typical seeing of about 0.9 arcsec in K-band. 

Four fields were observed for A2597 and three fields for 
S159. These fields were selected to lie within areas known 
to have extended Ha emission (J05). The observed spectral 
and spatial resolution, as measured from telluric lines and 
standard star observations, is summarised in Table (2] 

2.1.1 Initial Reduction and Slit Definition 

The reduction of the data was done using a combination 
of the ESO SINFONI pipeline recipes (SINFONI pipeline 



version 1.7.1 and GPL version 3.6.1), ECLIPSE (|Devillardl 
[200 r) and a set of dedicated IDL procedures. From the SIN- 
FONI pipeline we obtain a masterdark frame, masterfiat 
frame, hot pixel map and a slit curvature model. Wavelength 
calibration, hot pixel removal, slit edge detection, distortion 
correction, sky removal and illumination correction as given 
by the pipeline were found to be inadequate for our pur- 
poses and therefore an additional set of reduction tasks was 
written in IDL. 

The reduction was carried out as follows. Source and 
sky frames are corrected for dark current and fiat fielded 
using the masterfiat and masterdark frames from the SIN- 
FONI pipeline. Having estimated the slit edges (by eye) the 
different slits are defined and cut out of the science frames. 
Each of these slits is then treated independently in the sub- 
sequent reduction steps. CCD artefacts are removed from 
the data. Hot pixels and those affected by cosmic rays are 
interpolated over. 

We correct for the spatial curvature of the slit optics 
as imaged on the detector by applying the curvature model 
obtained by the SINFONI pipeline using the ECLIPSE task 
warping (Dcvillard 2001). Slit columns which do not contain 
information across the full wavelength range are removed. 
This also removes the overlap between different slits as im- 
aged on the CCD. 

The spectra are wavelength calibrated using a set of 19 
identified telluric OH lines in t he wavelength range 1.95- 
2.30 jUm (|Rousselot et all I2OO0I ). The output wavelength 
scale is set to 2.45 x 10~* /im per pixel thereby Nyquist sam- 
pling the data. 

2.1.2 Sky Removal 

The K-band night sky is variable on short time scales. We 
have rather long exposure times, as compared to the vari- 
ations in the sky background. This means that there is a 
complicated relationship between the sky contribution in 
our source frame and the sky observed in our corresponding 
sky frame. This is readily observed by subtracting two sky 
frames taken directly after each other, and leads to system- 
atic residuals of up to 10% in the peak heights of telluric 
lines. A scaling between the source and its corresponding 
sky frame thus needs to be performed in order to decrease 
these residuals. The standard sky scaling performed by the 
ESO SINFONI pipeline reduces these residuals to about 5% 
and the special SINFONI pipeline sky correction utility re- 
duces the residuals further to about 4%. In both cases it was 
noted that the sky removal suffered from poor wavelength 
matching between the sky and source frames due to flexure 
of the instrument. 

In our approach we remove the sky emission after de- 
tailed wavelength calibration using the telluric lines. We a 
apply a simple scaling function to the sky frame before sub- 
tracting it from the source frame. This scaling function con- 
sists of a constant and a small, linear, wavelength dependent 
factor. The constant is determined from the relative heights, 
above the continuum, of the telluric lines and assumed to 
hold at 2.1 fim. The small, linear, wavelength dependent 
factor is the slope of a linear fit to the ratio of the object 
spectrum and its corresponding sky spectrum. 

The full wavelength dependent behaviour of the sky 
emission between an object and a sky frame is often more 
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complex than the simple linear function used. Here we are 
only interested in line emission and as such a small resid- 
ual gradient in the continuum emission does not affect the 
analysis performed below. Our method leads to residuals 
that are ^2% in the peak heights of telluric lines. This is a 
significant improvement over the other methods mentioned 
above. In the final analysis of the line emission we checked 
our results carefully for telluric line residuals and removed 
wavelength regions affected by these from our analysis. 

2.1.3 Illumination Correction 

After correcting for any residual distortion we collapse the 
sky and the sky subtracted source spectra into cubes with 
pixel size (0.125" ,0.125" ,0.000245 ^m). It is known that 
SINFONI, after all reduction steps described above, still has 
a varying illumination across its FOV and that this illumina- 
tion is a function of wavelength (J. Reunanen priv. comm.). 
This is mostly due to a difference in the illumination of the 
various slitlets and most easily observed in the sky cubes. 
Defining a reference slitlet in the sky cubes we determine 
the variation in illumination across the FOV for each wave- 
length. We then correct for this variation in the sky sub- 
tracted source cubes. The correction is typically less than 
10% and particularly affects wavelengths below 2.1 /im. 

2.1.4 Flux Calibration 

Flux calibration is carried out using one or, if available, mul- 
tiple standard star observations per night. The standard star 
observations are reduced in the same way as outlined for the 
object observations above. All standard stars observed were 
either O or B stars, and brighter than 8th magnitude in 
the K-band. The atmospheric transmission function was de- 
termined by dividing the spatially summed standard star 
spectrum with a black body spectrum of the appropriate 
temperature. The absolute flux scale is set by using 2MASS 
K-band magnitudes, these are accurate to 0.05 magnitude 
in K-band. 

2.1.5 Extracting the line emission 

Following flux calibration the source cubes are combined. 
The northern and southern edges of the exposures for the 
different fields overlap well. The eastern and western edges 
overlap less well leading to a higher noise here. The most 
northern and southern slitlet have very low signal to noise 
and were removed from the data. Any remaining telluric 
emission is removed by defining off-source regions. These 
are marked by the dotted lines in Figs. [1] and O 

Continuum emission is removed by fitting the contin- 
uum in the immediate neighbourhood of the science line. 
The continuum subtracted line feature is fit by a single 
Gaus sian function, using the mpfitpeak routine (jMarkwardtl 
I2OO9I ) within IDL. It is observed that a single Gaussian al- 
ways provides a good description of the observed line pro- 
file. The line flux, centre and width are determined from this 
Gaussian fit. For selected regions line profiles and Gaussian 
fits to them are shown in Appendix [C] The errors quoted 
for the fitted line properties are based on Monte-Carlo sim- 
ulations. 



2.1.6 Constructing the line maps 

In order to visualise the surface brightness of the line 
emission we performed a Gaussian smoothing of four 
pixels full width at half maximum (FWHM) in both 
the spatial (4 pixels=0.5") and the spectral plane (4 
pixels=9.80x 10~* ^m). To visualise the kinematics of the 
line emitting gas a two pixel FWHM Gaussian spatial 
smoothing and no spectral smoothing was found to be ade- 
quate for A2597, whereas for S159 a two pixel FWHM Gaus- 
sian smoothing in both the spatial and spectral planes was 
performed. The degradation of the spatial and spectral res- 
olution as a function of the smoothing kernel used is given 
in Table [2I The corresponding noise is given in Tables [3] and 

H 

Surface brightness maps for all lines that could be 
mapped on a pixel to pixel basis are shown in Appendix 
1X1 and [B] For A2597 the northern field is not shown as 
the signal to noise here is inadequate to show the emission 
on the same spatial resolution as the central and southern 
fields. Velocity and velocity dispersion are shown only for the 
strongest detected ionised and molecular gas line. We note 
that the velocity structure observed in all detected emission 
lines agrees with that shown for these lines. 

2.2 X-ray Data 

Cool core clusters were first discovered by analysing their 
X-ray emission. These observations lead to the still un- 
resolved cooling flow pr oblem for the hot X-ray gas (e.g. 
[Peterson fc FabianI [2OOT ) . In this paper we are concerned 
with the cooler HII and H2 gas phases and will not focus on 
the cooling flow problem related to the hot X-ray gas. How- 
ever, in Section 8 we will investigate whether there is a spa- 
tial correlation between the X-ray emitting gas and cooler 
gas phases. In order to do so we have obtained all available 
X-ray data from the CHANDRA archive. The A2597 image. 
Fig. [1] is a co-add of three separate observations having a 
combined exposure time of 153.7 ks (project codes 7329; 
6934; 922). The S159 image. Fig. [2] consists of only one 
shallow 10.1 ks observation (project code 1668). 

CHANDRA data for A2597 and S15 9 has previously 
been published bv lMcNamara et al.l (|200ll ) and J05. A very 
notable difference in the X-ray emission for the two clusters 
is that the peak emission in A2597 is well aligned with the 
nucleus of the BCG, whereas in S159 the peak emission is 
about 8" north of the BCG nucleus. 

2.3 Radio Data 

Out of the many heat sources proposed, AGN feedback has 
received the most attention in recent years. The observed 
anti-correlation between X-ray and radio emission, also re- 
ferred to as X-ray cavities and Radio bubbles, has led to 
models in which the AGN outflows interact strongly w ith its 
surrounding medium (e.g. [Sutherland fc Bickneiill2007h . The 
kinetic energy carried by these outflows has been calculated 
based on these X-ray cavities and recent results show that 
the mechanical power of the jet that created the X-ray cavi- 
ties can be orders of magnitude large r than its radio inferred 
radiative power (jBirzan et al.l |2004| . l2008l : |d unn fc FabiarJ 
l2006h . In Section 8 we will compare our SINFONI results 
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with high resolution radio maps to investigate how the cur- 
rent AGN outflows interact with the cooler gas in the cores 
of A2597 and S159. 

A2597: A V LA 8.4 GHz map of A2597 was obtained 
from C. Sarazin jSarazin et al.|[l995l '). This map has a beam 
size of 0.26" X 0.21". The one sigma noise is 50 /xjy beam"^. 
It was noted that there is a significant offset of ~0.1 sec- 
onds in Right Ascension as compared to the 2MASS posi- 
tion of the BCG. Two more 8.4 GHz maps ha ve been pub- 
lished (jBirzan et al.ll2008l : lOonahue et aLlbood ). These have 
a much better agreement with the 2MASS position. We thus 
conclude that this offset is an error and shift the 8.4 GHz 
map accordingly. Detail ed radio maps at lo wer frequencies 
have been pubUshed bv lClarke et~al] (|2005l ) and show that 
there is more low level radio emission present than appar- 
ent from the 8.4 GHz map. However, the 8.4 GHz map does 
give a good indication of the current AGN outflows. A much 
higher resolution radio map at 1.3 and 5.0 GHz using very 
long bas eline array (VLBA) interferometry has been pub- 
lished bv lTavlor et al.l (|l999l ). These observations show that 
the current jet has a position angle (PA) of 70 degrees. 

S159: Archival VLA 8.4 GHz observations of S159 
(project code: AB1190) were reduced with the NRAO 
Astronomical Image Processing System (AIPS). The B- 
configuration observations were taken in single channel con- 
tinuum mode with two IFs centred at 8435 and 8485 MHz. 
The total on source time was 103 min. The data was fiux cali- 
brated using the primary calibrator 013 7-1-331. We u s ed the 
Perley & Taylor 1999 extension to the iBaars et all (Il977l ') 
scale to set the absolute flux scale. Amplitude and phase 
variations were tracked using the secondary calibrator 2314- 
449 and applied these to the data. The data was imaged 
using robust weighting set to 0.5, giving a beam size of 
3.26" X 0.67". The one sigma map noise is 25 ^tjy beam~^. 
Radio maps of S159 at 1.4 , 5.0 an d 8.4 GHz were previously 
published bv lBirzan et all (120081 ). 



3 ABELL 2597 - GAS DISTRIBUTION 

Four 8"x8" fields were observed on and surrounding the 
BCG PGC071390 in A2597, see Fig.ffl The integration time 
for each exposure is 600 seconds. The central field, which 
includes the nucleus of PGC071390, contains 13 exposures. 
The south-eastern (SE) and south-western (SW) field con- 
tain 8 and 15 exposures respectively. The northern field con- 
tains 13 exposures. The overlap region between the central 
and southern fields is sufficient for the line emission to be 
mapped without problems. However, the SE and SW fields 
do not completely overlap everywhere. In various locations 
along the overlap area there are small gaps between the fields 
of one to two pixels (1 spatial pixel=0.125"). We interpo- 
lated these before mapping the emission. Despite this, due 
to the increased noise at the east, west edges of each field, 
this overlap region (about 10 pixels in width) between the 
southern fields has a rather poor signal to noise. The north- 
ern field is offset by about 6" from the central field. 

A four pixel spatial and spectral smoothing was applied 
to the data prior to fitting the lines. A single Gaussian func- 
tion provides a good fit to the observed line profile. Surface 
brightness maps for all other lines that could be mapped on 
a pixel to pixel basis are shown in Appendix |X] The north- 



ern field is not shown in these images because the signal to 
noise here is inadequate to show the emission on the same 
spatial resolution as the central and southern fields. 

3.1 Molecular gas 

The integrated line fluxes for all lines detected within the 
observed fields for A2597 are given in Table [S] AU H2 1-0 
and H2 2-1 S-transitions redshifted into the K-band (1.95- 
2.40 lira) are detected. A fiux value for the H2 2-1 S(4) line 
has been omitted due to uncertain continuum subtraction. 
The H2 2-1 S(5) and the Br 5 line are too closein wavelength 
to be separated by our observations. None of the H2 3-2 S- 
transitions were detected. As an example of the fidelity of 
the data we show the full K-band spectrum of the nuclear 
region and the south eastern filament in Fig. |4] 

The H2 surface brightness maps all show the same struc- 
ture. As an example of the molecular gas distribution we 
show the surface brightness map for the H2 1-0 S(3) line 
in Fig. |S] This map clearly shows that the peak of the 
molecular gas emission coincides with the stellar nucleus of 
PGC071390. Two extended gas structures away from the 
nucleus are observed. One extends north from the nucleus 
and hence we will refer to this structure as the northern fil- 
ament. The second structure extends from the north-east to 
the south-west in the SE field, just south from the nucleus 
and hence we will refer to this structure as the southern 
filament. 

We observe that the surface brightness of the molec- 
ular gas varies rather smoothly within the nuclear region. 
However, from higher spatial resolution HST imaging by 
iDonahue et~al] (|2000l . hereafter DOO) we know, that the 
molecular and ionised line emission in the central 4"x4" 
is concentrated in narrow clumpy, filamentary structures. 
Here we do not have the resolution to resolve these struc- 
tures. We do note some enhanced emission features, embed- 
ded within the central field, extending to the north and east 
away from the nucleus which are roughly coincidental with 
some of structures observed by DOO. 

The northern filament extends at least up to the north- 
ern edge of the central field, i.e., 6" (9 kpc) from the nucleus. 
This is well beyond the region in which molecular emission 
was detected by DOO. Using deep K-band longslit spectra 
J05 have previously observed that the H2 emission extends 
at least up to 20 kpc towards the north from the nucleus. 
We will see below that molecular gas can still be found in 
the northern field observed by us, i.e., at a distance of about 
22 kpc from the nucleus, thus confirming the J05 results. 

The southern filament is clearly detected in the emission 
of the stronger lines. This southern filament has not been 
observed in DOO, but J05 also find molecular gas south of 
the nucleus (see their figures 8 and 11). The extent of the 
northern and southern filaments observed here is bounded 
by the edges of the observed fields, and it is likely that these 
continue beyond the regions mapped by us. 



3.2 Ionised gas 

The Pa ct line is redshifted into the K-band for both galaxy 
clusters studied here. The line is redshifted into a region 
of poor atmospheric transmission, but it is the strongest 
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ionised gas line by far in our spectra and unambiguously de- 
tected in both clusters. In A2597 the Pa a emission globally 
follows the H2 emission closely, Fig.|6l Within the nuclear re- 
gion enhanced emission is again observed towards the north 
and east. These features are roughly coincidental with the 
emission line filaments observed in DOO. Beyond the nuclear 
region the emission extends along the northern and southern 
filaments, peaking in the same locations as the H2 emission. 

We also detect the Br 7, Br S and Fe II (1.8100 fim 
rest wavelength) lines. The Br 5 line is blended with H2 2-1 
S(5) and these cannot be disentangled directly by our ob- 
servations. In the central field we find that Br 7/Pa a — 
0.082. This ratio agrees with the dust-free Case B recombi- 
nation ratio of these lines for ne = lO'^ cm~^ and T = lO** K 
llOsterbrock fc Ferland|[2006h . The Case B scenario then im- 
plies that Br 7/Br 5 — 1.5, and we use this ratio to disen- 
tangle the Br 5, H2 2-1 S(5) blend. In the nuclesar region 
we find that Br 5 and H2 1-0 S(5) are of similar strength. 

A small dust la ne has been observed i n the nuclear re- 
gion of A2597 (DOO: lKoekemoer et al.lll999l ). We have inves- 
tigated whether differential extinction in the K-band may 
affect our emission line ratios. From the above value for the 
Br 7/Pa a ratio we find that differential extinction is unim- 
portant in the K-band. This is co nfirmed by d eep optical 
spectroscopy of A2597 by VD97 and lBaked l|2005l ). They find 
a V-band extinction Av ~ 1 across the nebulosity. Assum- 
ing standard galactic dust (Rv=3.1) an Ay ~ 1 translates 
in to Ak ~ 0.1. This amount of extinction is negligible. 

The Fe II (1.81 /xm) line is redshifted to the short wave- 
length edge of our observed spectrum. It is unambiguously 
detected in the nuclear region. The decrease in the Fe II 
emission outwards from the nucleus, in both intensity and 
dispersion, appears to be much faster than for either the HII 
lines or the II2 lines. The HII emission drops by a factor of 3 
and the H2 fiux drops by a factor 4 from the nuclear region to 
a region just north of the nucleus. The Fe II emission drops 
by a factor of 10 for the same regions. If the Fe II emission 
has a different origin than hydrogen lines, for example if it 
is preferably coming from the AGN and the associated jet 
instead of the gaseous filaments, this may explain the differ- 
ence. Our observations do not have the spatial resolution to 
investigate this in detail. 

We have searched our spectra for the presence of even 
higher ionisation lines, such as the Si VI (1.9634 fim) line, 
which one would expect from typical hard AGN spectra. 
None of these higher ionisation lines were detected. This 
once more confirms the LINER nature of PGC071390. It 
may also indicate that the active nucleus is not the main 
source of ionisation of the gas observed in the core of A2597. 
Alternatively it would have to have an atypically soft ionis- 
ing spectrum. 



4 ABELL 2597 - GAS KINEMATICS 

A single Gaussian function gives a good description of the 
observed line profiles, see Appendix [C] From the Gaussian 
fits of these line profiles we obtain information about the 
kinematical structure of the molecular and ionised gas in 
A2597. The velocity, with respect to the systemic velocity 
of PGC071390, and the velocity dispersion of the gas have 
been derived for all emission lines. These all show the same 



global kinematical structure. The velocity and velocity dis- 
persion maps shown below differ from the surface brightness 
maps in that only a two pixel spatial smoothing has been 
applied, as supposed to a four pixel spatial and spectral 
smoothing. This is done to preserve as much of the velocity 
structure as possible and provides us with a velocity resolu- 
tion of 38 km s~^. 



4.1 Molecular gas 

The molecular gas in A2597 shows a wealth of small scale 
kinematical structure. Velocity and velocity dispersion maps 
were made for all H2 lines. All show the same kinematical 
structure on all scales observed. As an example of this struc- 
ture we have displayed the velocity and velocity dispersion 
maps for the H2 1-0 S(3) line in Figs. [7] and [S] The nuclear 
region contains a blueshifted and a redshifted gas compo- 
nent at about ±80 km s~^. This is seen more clearly if we 
place a pseudo slit with a width of 1" and a PA of 105.5 
degree, centred 1 kpc south of the stellar nucleus. The cor- 
responding position-velocity diagram is shown in Fig. 1151 
The velocity structure observed in Figs. [7] and [T^] is reminis- 
cent of gas rotating around the nucleus and does not appear 
to be related to an expanding shell or AGN outflows. 

The average velocity of gas in the nuclear region is ap- 
proximately zero with respect to the systemic velocity of 
PGC071390 (z=0.0821, VD97). This shows that the gas here 
is situated at or near the stellar nucleus. The reason for plac- 
ing the pseudo slit slightly south of the nucleus is because 
east of the nucleus there is a small strongly redshifted feature 
at -1-150 km s~^. Whether this feature is part of the global 
gas flow or a single event is unclear. It shows up prominently 
in all velocity maps. Projected on to the sky, the feature 
appears to be coincident with the north-eastern radio jet of 
PKS2322-12 the radio source in PGC071390, see Fig. [7] The 
filamentary structures extending towards the north and the 
south from the nucleus show smooth velocity gradients and 
these will be discussed in more detail below. 

The velocity dispersion of the molecular gas also shows 
interesting structure. Globally the dispersion of the gas de- 
creases with distance from the nucleus. It drops from an 
average of about 220 km s~^ in the nuclear region to about 
100 km s~^ a few kpc north and the south of the nucleus. 
The velocity dispersion is very high in two narrow struc- 
tures extending towards the east and south of the nucleus. 
The two-dimensional data allows us to determine the area 
which is disturbed to be an elongated structure of about 
2 kpc by 5 kpc oriented at a PA of about 45 degrees. 

Projected on to the sky these high dispersion structures 
appear to run along the inner. South East edge of the curved 
radio lobes of PKS2322-12, see Fig. [S] If we interpret the 
lobe morphology as a Wide Angle Tail, caused by the rela- 
tive motion of the AGN and the external medium, then the 
dispersion map illustrates for the first time the turbulent 
wake expected from this motion. Alternatively, the region of 
maximum dispersion, at PA~70 degrees from the nucleus, 
may represent the interaction of the current, VLBI radio jet 
wit h the surrounding med ium, as has been seen in Centaurus 
A l|Neumaver et al.ll2007l ). In this picture we must assume 
that the counter-jet has been defiected near the nucleus to 
the South, causing the high dispersion region and radio lobe 
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in this direction. There is, however, no evidence for a major 
kinematic disturbance at the point of deflection. 

The highest velocity dispersion is found for the small, 
strongly red-shifted feature east of the nucleus. This high ve- 
locity and dispersion for this feature can be explained if this 
is gas that is entrained within the AGN outflow. The fea- 
ture alig ns well with the cu rrent, projected jet axis (PA=70 
degrees, iTavlor et alll 19991 ) . 

4.2 Ionised gas 

Velocity and the velocity dispersion maps for the ionised gas 
in A2597, as traced by the Pa a line, are shown in Figs. [7] 
and [S] We observe two key features when we compare the 
Pa a and H2 derived kinematics. Firstly, globally we find 
that the Pa a derived kinematics follows the H2 derived 
kinematics tightly. Secondly, the velocity dispersion of the 
Pa a emitting gas appears on average to be slightly higher 
than the H[2 emitting gas, especially in the nuclear region. 

It may be possible that on scales below the resolution 
of our observations the ionised gas has a different distribu- 
tion than the H2 emitting gas. This may be especially true 
in the nuclear region where the active nucleus appears to 
be strongly interacting with the gas. The position-velocity 
diagram shown in Fig. [15] also shows that the ionised g 
traced by the Pa a and Fe II lines, reaches slightly higher 
velocities in the nuclear region. DOO show that within the 
nuclear region the H2 and HII gas has a very complicated 
and disturbed morphology and it is difficult to say how well 
these two trace each other on small scales here. 

The kinematics of the molecular and ionised gas for 
A2597 derived here agrees well with p revious lon g slit in- 
vestig ations by JOS and Heckman et al.l |l98 9). Q'Dea et al.l 
l| 19941 . hereafter 094) detected HI in absorption against the 
radio continuum source PKS2322-12 in A2597. The absorp- 
tion observation represents a line of sight of a few arcsec 
along the central radio source. They find a spatially resolved 
broad HI component with a ~174 km s~^ and a narrow un- 
resolved HI component with a ~93 km s"'^ at the position 
of the nucleus. The width of the broad component is some- 
what smaller than the width observed in HII and H2. 094 
find that the widths are consistent if one takes into account 
that the HI absorption measurements only sample the gas 
in front of the radio source, whereas the HII and H2 mea- 
surements sample all of the gas along the line of sight. 

As in our data 094 find a narow and a broad com- 
ponent, but the relative strength of narrow component is 
much stronger in their observation. We do not see the nar- 
row component on the nucleus. The dominance of the narrow 
component in the HI observations is probably caused by the 
1/Tb dependence of the HI absorption, as pointed out by 
094. Tb being the spin temperature of the HI gas. In the HI 
absorption spectra the cold gas at large radii in front of the 
nucleus is probably over-represented relative to the HII and 
H2 emission spectra. We conclude, as do 094, that there is 
no evidence for a kinematically distinct HI component. 

4.3 Filaments 

In Fig. [13] we show the surface brightness, velocity and ve- 
locity dispersion along the two filamentary structures we 



identified in our observations of A2597. The regions used for 
this investigation are marked by the green and red squares 
in Fig. O The black points in Fig. [T3] correspond to green 
squares and the red points to the red squares. Following 
the northern filament from slightly south of the nucleus to- 
wards the northern edge of the central field we find that 
the Pa Q/H2 1-0 S(3) is approximately equal to 0.75 in the 
nuclear region and rapidly increases to unity outwards. The 
northern filament shows a smooth velocity gradient from 
south to north across the nucleus, as the velocity decreases 
from -1-50 km s~^ to -50 km s~^ 2" north of the nucleus. 
At this point the velocity gradient reverses and the velocity 
increases again to -)-50 km s"'^ towards the northern edge of 
the central field. 

Velocity gradients and even reversals for this filament 
may be explained in terms of bending and stretching of the 
filament, perhaps due to a combination of its proper mo- 
tion and gravitational forces. However, it is more likely that 
we observe multiple filaments, each with its own character- 
istic motion, along the line of sight. Our data shows that 
the eastern part of the northern filament is predominantly 
blue shifted whereas the western part is red sh ifted. Higher 
spatial resolution images taken with HST by lO'Dea et al.l 
(2004) and Oonk et al. (in prep.) show evidence that the 
northern filament observed here consists of at least two fil- 
amentary structures. We thus favour the latter explanation 
for the observed velocity structure of the northern filament 

This interpretation also agrees with what is observed for 
more nearby galaxy clusters such as Perseus and Centaurus, 
where a multitud e of long, th in filaments ar e observed along 
the line of sight (iFabian et al. 2008; Crawford et all I2OO5I : 
iHatch et al.ll2005l ). The narrow spatial and velocity range 
observed here for the filaments however still suggest that 
any substructure in it will likely have a common origin. If 
the gas observed in the northern filament is connected to the 
gas detected in the northern field its velocity continues to 
increase to about -f 150 km s"'^, as also shown by J05. From 
the J05 observations it appears that the gas in the central 
field is joined smoothly with that in the northern field, in 
terms of both surface brightness and dynamics. 

The velocity dispersion along the northern filament de- 
creases smoothly from 220 km s~^ to 100 km s"'^, from the 
nucleus to the edges of the central field. This decrease is 
fastest near the nucleus and slows down beyond 3 kpc north 
of the nucleus. This point may mark a change in the influ- 
ence of the AGN upon the dynamical state of the gas. 

The southern filament has a much lower surface bright- 
ness and is hence detected at a lower signal to noise. Varia- 
tions along this filament are thus more difficult to detect. 
Following this filament from the north-east (NE) to the 
south-west (SW) we find that the surface brightness is high- 
est at its NE edge whereafter it decreases slightly and be- 
comes approximately constant. The velocity decreases from 
-1-50 km s~^ to about -40 km s~^. The velocity dispersion re- 
mains constant at about 100 km s~^ along the filament. We 
will discuss the stability of the observed filaments in more 
detail below. 
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5 SERSIC 159-03 - GAS DISTRIBUTION 

Three 8" x 8" fields were observed on and surrounding the 
BCG ESO291-G009 in S159, Fig. [3 The integration time for 
each exposure is 600s. The south-eastern (SE) and south- 
western (SW) fields contain 8 and 9 exposures respectively. 
The northern field contains 8 exposures. The SE field con- 
tains the nucleus of ESO 291-G009. There is no overlap be- 
tween the three fields observed. A four pixel spatial and 
spectral smoothing was applied to the data prior to fitting 
the lines. A single Gaussian function provides a good fit to 
the observed line profiles. Surface brightness maps for all 
detected emission lines that could be mapped on pixel to 
pixel basis are shown in Appendix iBl 



5.1 Molecular gas 

Only two out of five H2 1-0 S-transition lines redshifted into 
the K-band are unambiguously detected for S159. These are 
the H2 1-0 S(l) and H2 1-0 S(3) lines. Their resulting inte- 
grated line fluxes are given in Table |6l The presence of the 
H2 1-0 S(2) line can also be confirmed, see also Fig. [S] This 
line cannot be fitted reliably due to a strong telluric line 
residual in the red wing of the line profile and as such a flux 
value has been omitted. Tentative evidence is found for the 
presence of the H2 1-0 S(0) and H2 1-0 S(4) lines. However, 
both lines lie in spectral regions of poor atmospheric trans- 
mission and are close to strong telluric features, preventing 
us from claiming detections. We have searched for the H2 2- 
1 S-transitions, but these and higher H2 transitions remain 
undetected in our spectra. Except for the above mentioned 
two H2 lines only the Pa a line is detected. Full K-band 
spectra of the nuclear region and the western filament are 
presented in Fig. O 

The two H2 surface brightness maps show the same 
structure and as an example of this we show the map for 
the H2 1-0 S(3) line in Fig. (9] This map clearly shows that 
the peak of the H2 emission is well aligned with the stellar 
nucleus of ESO0291-G009. A filament of molecular emission 
extends north-east from the nucleus up to the edge of the 
SE field. We will refer to this structure as the northern fil- 
ament. Clumpy emission extends towards the west and the 
south of the nucleus, up to the edges of the SE field. 

The SW field shows a strong filament of gas hav- 
ing an east-west elong ation, originally discovered by 
ICrawford fc FabianI J 19921 ). We will refer to this structure as 
the western filament. The northern field shows two features 
(i) low signal to noise clumpy emission in the southern and 
central part of the field, and (ii) a stronger, somewhat larger 
emission feature at the northern edge of the field. Both fea- 
tures are treated in more detail below. Whether the clumpy 
emission observed in the northern field is part of the north- 
ern filament observed in the SE field is not clear from our 
observations. 

The spatial extent of the gas for both the western and 
northern filament is bounded by the edges of our observed 
fields and it is likely that these filaments continue beyond the 
regions mapped here, as seems to b e the case from narro w- 
band Hq-|-[NII] imaging by J05 and [Hansen et all l|2000h . 



5.2 Ionised gas 

Pa a is the only ionised gas line detected in our K-band 
spectra for S159. It is redshifted into a region of poor atmo- 
spheric transmission. Strong Pa a emission is found along 
the northern and western filaments. Fig. [9] We detect Pa a 
in all places where H2 emission is detected. A noticeable dif- 
ference concerns the nuclear region. Almost no Pa a emission 
appears to be associated with the stellar nucleus. As we will 
see below, some Pa a emission is found here, but there is 
a strong decrease of it relative to molecular hydrogen emis- 
sion. The Pa a/H.2 1-0 S(3) ratio is observed to drop from 
about 0.7 in the filaments to about 0.2 in the nuclear region. 
We note that the nuclear region has a rather high spectral 
noise, due to the strong stellar continuum, which affects our 
ability to detect emission lines here. The detection of the 
Pa a line is furthermore complicated by it being in a region 
of poor atmospheric transmission. 

Pa a is also present within the northern field and, like 
the H2 gas, appears to be clumpy. The strong emission fea- 
ture observed in H2 towards the northern edge of this field 
is also confirmed by Pa a emission. 

In SI 59 the Pa a line is the only ionised gas line de- 
tected. An estimate of differential extinction within the K- 
band can th us not be made directly from our observations. 
I Baker! ()2005l ) finds finds little variation from Ay ~ 1 across 
the nebulosity. Applying the same arguments to S159 as for 
A2597, we conclude that differential extinction in the K- 
band is negligible for this cluster. 



6 SERSIC 159-03 - GAS KINEMATICS 

A single Gaussian function gives a good description of the 
observed line profiles, see Appendix [C] The velocity, with 
respect to the systemic velocity of ESO 0291-G009, and the 
velocity dispersion of the gas have been derived for all emis- 
sion lines. These all show the same global kinematical struc- 
ture. The velocity and velocity dispersion maps shown below 
differ from the surface brightness maps in that only a two 
pixel spatial and spectral smoothing has been applied. This 
is done to preserve as much of the velocity structure as pos- 
sible and provides us with a velocity resolution of 51 km s~^. 

6.1 Molecular gas 

The kinematical structure observed in the H2 1-0 S(l) and 
H2 1-0 S(3) maps is the same. As an example of this struc- 
ture we show the velocity and velocity dispersion of the H2 
gas as traced by the H2 1-0 S(3) fine in Figs. [10] and [TT] The 
nuclear region contains a blue- and redshifted gas compo- 
nent at about ±120 km s~^. This velocity structure is rem- 
iniscent of gas rotating around the nucleus. However, the 
two gas filaments coming in from the north-east and west 
towards the nucleus may also explain the observed velocity 
structure. The velocity along the filaments will be treated 
in more detail below. The gas extending towards the west 
and south from the nucleus appears to be predominantly 
blueshifted. The average velocity of the gas in the nuclear 
region is e qual to the syste mic velocity of ESO 291-G009 
(2=0.0564, iMaia et al.|[l987l ) showing that the gas is situ- 
ated at or near the stellar nucleus. 
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Globally the dispersion of the gas is low and decreases 
with distance from the nucleus. The dispersion of the gas in 
the nuclear region is about 230 km s~^, but drops rapidly to 
about 90 km s~^ along the two filaments. This is similar to 
what is observed in A2597. In projection the high dispersion 
structure around the nucleus appears coincidental with the 
radio jets of ESO0291-G009, see Fig.[TTl The increase in the 
velocity dispersion here again indicates that the radio jets 
are stirring up the gas. 

From Figs. llOland llll it is difficult to draw conclusions 
on the average velocity and velocity dispersion of the clumpy 
low signal to noise emission in the northern field. From the 
previous investigation by JOS it appears that the gas ob- 
served here is connected to the strong filament extending 
north from the nucleus. Below, we will see that the clumpy 
emission has velocities varying between -20 and -60 km s~^ 
and a velocity dispersion less than 100 km s~^. If this emis- 
sion is connected to the filament extending north this implies 
that the line emission continues to decrease with distance to 
nucleus as was also shown by J05. 



6.2 Ionised gas 

Velocity and velocity dispersion maps for the ionised gas in 
SI 59 as traced by the Pa a line are shown in Figs. [10] and 
111! As for A2597 we find that the Pa a derived kinematics 
follows that of the H2 lines closely. The only exception in 
S159 being the high dispersion nuclear region observed in H2 
emission, which appears to be missed by the Pa a emission. 
The high noise in the nuclear region of S159 combined with 
the poor atmospheric transmission in the wavelength range 
of the Pa a line makes the fit of this line here rather difficult. 



6.3 Filaments 

In Fig.[T4]we show the surface brightness, velocity and veloc- 
ity dispersion along the two filaments observed in S159. The 
regions used for this investigation are marked by the green 
squares in Fig. [S] Following the northern filament from the 
north-eastern edge in the SE field southwards toward the nu- 
cleus and subsequently to the eastern edge of the SE field, 
we find that the Pa a/H2 1-0 S(3) decreases strongly in the 
nuclear region and is approximately constant outside of it. 
The filament has a smooth velocity gradient. North-east of 
the nucleus the velocity decreases very slightly with distance 
from about -1-130 km s~^ to about -1-90 km s~^. Across the 
nucleus the velocity changes from about +100 km s~^ to 
-100 km s~^. Whether this velocity reversal is due to rota- 
tion or due to the filamentary structure of the gas can not 
be clearly distinguished. 

The velocity dispersion of the gas in the northern fil- 
ament is low everywhere, except within the nuclear region. 
All detected lines show an increase in the velocity dispersion 
near the nucleus, but the increase in the H2 lines is much 
higher than for the Pa a line. The decrease in the dispersion 
east of the nucleus is difficult to measure due to low signal 
to noise here. 

The surface brightness of the western filament has two 
prominent peaks about 11" from the nucleus. It has a very 
smooth velocity gradient along the filament, from about - 
100 km s~^ to about +20 km s~^. This velocity structure 



agrees with the possibility that it is connected to the nu- 
clear region. The dispersion of the gas in this filament is low 
everywhere. All three lines detected show the same flux and 
velocity structure along the western filament. 



7 PHYSICAL CONDITIONS IN THE WARM 
MOLECULAR GAS 

Molecular hydrogen emission can be excited by various 
processes, (1) thermal excitation produced by kinetic en- 
ergy inj ection into the gas due t o for example shock 
heating (|Hollenbach fc McKei 1 19891 '). (2) fluorescence by 
soft-UV photons, i.e. Photo dissociation regions (PDR ) 
(|Black fc van Dishoecklll987l : [Sternberg fc Dalgarnolll989l '). 
(3) high energy X-ra y photons, i.e. X-ray dissociation regions 
XDR) |Malonev et al. 1996 ) and (4) high energy particles 
Lepp fc McCravlll983l : iFerland et al.l 12009'). If the density 
of the molecular gas is above the critical density the gas is in 
local thermal equilibrium (LTE). It is then not possible to 
distinguish between the various H2 excitation mechanisms 
and one would observe a thermal H2 spectrum where the 
excitation temperature is equal to the kinetic temperature 
of this gas. For the photon and particle processes mentioned 
above the gas can be thermalised, via heating through sec- 
ondary (suprathermal) electrons. 

In order to investigate the II2 excitation mechanism we 
have constructed H2 excitation diagrams for a number of 
regions in A2597. Seven regions, A1-A7, were selected based 
on having distinct physical properties in terms of either their 
emissivity, spatial location or kinematical structure, see Fig. 
(3] Similarly seven regions, B1-B7, were selected in S159, see 
Fig. [3] Excitation diagrams were not made for S159 since we 
only have reliable detections for two H2 lines. 

All line detections for regions A1-A7 and B1-B7 are 
presented in Appendix[Cl The line proflles are well described 
by a single Gaussian function. For the lines detected in these 
regions their integrated fluxes are given in Tables Ffl and 1 101 
Kinematical information for these regions are presented for 
the H2 1-0 S(3) and Pa a lines only and these are given 
in Tables [8] and 1111 Gas temperatures and masses for the 
selected regions are given in Tables [9] and 1121 We will first 
shortly describe the selected regions below. 



7.1 A2597: Selected regions 

In the central field we have selected four regions, Al, A2, 
A3 and A6. Region Al corresponds to the nuclear region, 
A2 samples the region where the northern filament connects 
to the nucleus. A3, just north of A2, samples that part of 
the northern fllament that appears to not be influenced (di- 
rectly) by the nucleus anymore. Lastly, in the central field 
a clump of strong emission is noted on the western edge of 
the field, which we selected as region A6. There is tentative 
evidence for a narrow ridge connecting A6 to the nucleus, 
but we cannot confirm this with the present observations. 
The emission in region A6 itself is also uncertain, due to 
the increased noise at the edge of the field. It is only ob- 
served significantly in two H2 lines and therefore we have 
not constructed an excitation diagram for this region. 

In the south-eastern field we have selected one region, 
A4, which captures most of the emission in the southern 
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filament. This region was selected such to avoid the noisy 
overlap region between the SE and SW fields. The strongest 
emission lines showed weak evidence in their surface bright- 
ness maps that the south-eastern filament may extend across 
the overlap region into the SW field. Region A5 was selected 
to test this. Significant line emission is found in this region. 
The kinematical structure of the gas observed in A5 is sim- 
ilar to that measured in A4 and thus a connection between 
the two regions seems plausible. 

The northern field also contains significant emission for 
the strongest emission lines. A systematic search for line 
emission in the northern field was performed using various 
binnings. Region A7 was selected to show that emission does 
exist at a significant level in the northern field. JOS previ- 
ously showed that molecular and ionised gas existed out to 
20 kpc north from the nucleus using long slit spectra. The 
H2 1-0 S(l), 1-0 S(3) and Pa a lines are all detected at the 
3.0 sigma level, and the H2 1-0 S(5) line is observed at the 
2 sigma level. All four detected lines show the same velocity 
structure. 

The velocity dispersion of the gas in A7 is more difficult 
to constrain. Of the four lines the H2 1-0 S(3) line is observed 
at the highest significance and has the most reliable fit. This 
line has a dispersion of about 60 km s~^ using a spatial 
smoothing of two pixels. The other lines have higher fitted 
velocity dispersions ranging from 70 to 110 km s~^, however 
within their large errors (40-60 km s^^) they agree with the 
H2 1-0 S(3) result. 

The positive velocity and low velocity dispersion of the 
gas makes it plausible that the gas in A7 is connected with 
the northern filament. This interpretation agrees with the 
J05 results. We conclude that molecular and ionised gas is 
present at least up to region A7, i.e. 22.5 kpc from the nu- 
cleus, in good agreement with the J05 results. There is ten- 
tative evidence from regions investigated north and south 
of A7 within the northern field that line emission is present 
there as well. 

7.2 SI 59: Selected regions 

In the SE field we have selected four regions, Bl, B2, B3 
and B4. Region Bl corresponds to the nuclear region. B2 
samples the tail of the northern filament as it connects to 
the nucleus. Region B3, just north of B2, is selected to sam- 
ple the filament where it is no longer (directly) influenced 
by the nucleus. Lastly, B4 is selected to contain the clump 
of emission just to the south of the nucleus. Whether this 
clump is part of a filament or not our observations can not 
confirm. 

As discussed above, the nuclear region Bl is signifi- 
cantly noisier than surrounding regions. The spectral noise 
is higher here by a factor two to three. The H2 1-0 S(l) and 
1-0 S(3) lines are still easily detected here. The H2 1-0 S(3) 
line is considerably stronger than the H2 1-0 S(l) line in this 
region as compared to any of the other regions. The Pa a 
line is detected at a much lower significance and at a much 
lower flux level. In the SW field we have selected one region, 
B5, which captures most of the western filament. 

We have selected two regions in the the northern field, 
B6 and B7, to investigate the low level clumpy emission here. 
Region B7 was selected to capture the strong clump of emis- 
sion at the northern edge of this field. Region B6 was selected 



to investigate the remaining emission. The summed spec- 
tra clearly show that line emission is present in the north- 
ern field. We can thus conclude that molecular and ionised 
emission is present at least up to 18.0 kpc from the nucleus. 
Whether the gas in the northern field is directly connected 
to that in the SE field cannot be confirmed by our observa- 
tions, although it seems plausible from the JOS results. The 
velocity and low velocity dispersion of the gas are such that 
this gas can be connected smoothly to that in the northern 
filament. 

7.3 Thermal excitation of the molecular gas 

In the case of a gas in LTE, assuming a ortho:para abun- 
dance ratio of 3:1, there is a simple relation between the flux 
F and the temperature Tu corresponding to the energy of the 
upper state of a line, F ^ hu AN ~ hu qAexp(—T u /Te^c ), 
(e.g., JOS; IJaffe et al.l I2OOII : IWilman et all |2002| . l2005l ). 
Normalising the flux of each H2 line flux by the flux of the 
corresponding H2 1-0 S(l) line, we find; 



^ ' \ Fsiv,A,g, ) ^ 

= (^) X (ru,i - r,,si) (2) 

\ J cxc / 



Here F is the flux of the line, A its transition proba- 
bility, V its frequency and g the statistical weight of the 
transition. If the molecular gas is in LTE the H2 emission 
lines will lie on a straight line in a diagram of ln(F) vs. Tu. 
Tcxc, the reciprocal of the slope, will then be the kinetic 
temperature of this gas. We have investigated relation (1) 
for the all regions selected in Abell 2S97 in which we have 
detected at least 3 H2 lines. We show in Fig. [12] that the 
H2 lines detected in these regions are well fit by a thermal 
excitation model, with an average temperature of about 
2300 K. The derived excitation temperature for each region 
is given in Table [T] 

Besides the best-fitting LTE model for the H2 line 
fluxes we plot three additional H2 models. These models 
are shown for qualitative comparison purposes only, and 
are not tuned exactly to our physical conditions, see Fig. 
1121 The best-fltting LTE model is given by the black solid 
line. The red dotted line is a low-dens ity UV fluorescence 
model by iBlack fc van Dishoeckl l|l987D that does not in- 
clude collisions (their model 14; n=3xl0^ cm~^, a temper- 
ature r=100 K and a UV intensity /uv=10^ relative to 
/tot ) . The blue dotted line is a hi g h-den sity U V fluorescence 
model by ISternberg fc Dalgarnd ([198^) which does include 
collisions (their model 2D; n=lxlO® cm~^, a temperature 
r=1000 K and /uv=10^ relative to /tot). Lastly th e gree n 
dotted line is the cosmic ray model bv lFerland et al.l (|200i), 
which was developed for the gas filaments observed in the 
Perseus cluster. 

The Black & van Dishoeck low-density UV and the Fer- 
land cosmic ray models have distinct features that make 
them deviate from a thermal model. Low-density UV fluo- 
rescence models tend to boost the higher S-transitions (2- 
1, 3-2,...) relative to the 1-0 S-transitions as compared to a 
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LTE model. The Ferland model is observed to boost the even 
H2 1-0 S-transitions relative to uneven H2 1-0 S-transitions 
as compared to a LTE model. The high-density UV model by 
Sternberg & Dalgarno simply shows that at high densities, 
collisions within the gas will cause it to become thermalised 
and thus the line ratios also produce a straight line in our 
excitation diagrams. We thus conclude, qualitatively, that 
out of the four models investigated here that a LTE model 
provides the best description of the data. 

As the molecular line ratios in A2597 appear to be in 
LTE this implies that the densi ty of this gas is near its criti- 
cal density, nna.crit ~ 10® cm^ (|Shull fc Beckwithlll982l l and 
is dominated by coUisional excitation. Information regard- 
ing the source of excitation is thus not obtainable from this 
data set. 

There is a trend that on average we find higher LTE 
temperatures of the molecular gas in the filaments as com- 
pared to the nuclear region. We note though that within 
errors the temperatures agree for all regions, except for A6. 
Neither the H2 1-0 nor the H2 2-1 S-transitions lie exactly 
on a straight line. If we use only pairs of lines like the H2 1- 
S(l), 1-0 S(3) line pair or the H2 1-0 S(3), 1-0 S(5) line pair 
to determine an excitation temperature, assuming LTE, we 
find on average a temperature that is a few hundred degrees 
lower or higher than when all lines are used. Typically the 
first pair gives a lower temperature and the second pair a 
higher temperature. 

If the H2 gas observed here is at its critical density 
then there is a pressure imbalance between the molecular 
gas p{H2) ^ nT = ricrit Texc ~ 10^ and the ionised gas has 
p{HII) nT = 10^x10* = 10^. That the molecular and 
ionised gas are not in pressure equilibrium has previously 
been suggested by JOS using similar arguments. 

For Sersic 159-03 only two H2 lines were reliably de- 
tected in all regions. If we assume that the molecular gas 
here is in LTE, we can calculate an excitation temperature 
for this gas. The resulting excitation temperatures for the 
selected regions are given in Table [TOl Again we find an av- 
erage temperature of about 2300 K for the II2 gas in the 
filaments. In the nuclear region the strong increase of the 
II2 1-0 S(3)/H2 1-0 S(l) ratio leads to a very uncertain 
and physically unrealistic temperature above the dissocia- 
tion temperature for H2 . This either means that our line fits 
overestimate this ratio or that a different excitation mecha- 
nism is at work here. 

In this work we have only considered the K-band H2 
lines. These H2 lines do not show strong deviations from a 
single temperature thermal model with T ~ 2000-2500 K. 
From recent Spitzer spectroscopy (Jaffe & Bremer in prep.) 
we find that another H2 gas component exists in the nu- 
clear regions of A2597 and S159 with a LTE temperature of 
about 300 K. The situation in these clusters is thus similar to 
the situation in the Perseus cluster where multiple temper- 
ature LTE comp onents have been invoked t o explain the H2 
line r atios fe.g.. IWilman et al.|[20o3 . l2005l : [ Johnstone et al.l 
[2003). 

Whether the need for multiple LTE components to H2 
line ratios hints at a multiphase gas, a difference in optical 
depth, or a different excitation mechanism is unclear cur- 
rently. A more thorough modelling of the molecular gas, 
including all of the measured of the H2 lines in the infrared 



and mid-infrared, for A2597 and S159 will be presented by 
us in a future paper. 

7.4 Luminosity of the Warm Molecular Gas 

If we assume that the II2 gas is in LTE and can be described 
by a single kinetic temperature of about 2000 K, then the II2 
1-0 S(l) line luminosity represents about 10 percent of the 
total II2 luminosity. This is estimated using a list of 312 II2 
lines, corresponding to all II2 emission lines with an intensity 
greater than 1 percent of H2 1-0 S(l) line flux for TsC 4000K. 
From the total integrated II2 1-0 S(l) line fluxes in Tables 
[5] and [6] we thus find i(H2, A2597) = 1.1x10*^ erg s"^ and 
i(H2, S159) = 1.2x10*^ erg s'^ within the fields observed 
by us. 

This is significantly below the total H2 luminosity es- 
timated by J05, i.e. i(H2)~10''^~'"' erg s"^ J05 calculate 
the total II2 luminosity in the same manner as we do here. 
The difference follows from two simple arguments. Firstly 
J05 assumed that the H2 1-0 S(l) line luminosity represents 
only about 1 percent of the total II2 luminosity. This is true 
for low-density UV excitation models such as the Black & 
van Dishoeck models. However, w hen the gas is in LTE (s ee 
above) or in a XDR environment (|Draine fc Wood"3l 19921 ) a 
fraction of 10 percent is found. The JOS results thus need to 
be corrected down by a factor of 10 as consequence of this. 
The second argument that JOS made is that the H2 emis- 
sion coexists with the Ha emission. This argument follows 
from the HII/H2 ratio which is observed to remain rather 
constant over large areas (JOS and this work). JOS detect 
Ha over a much larger area than the fields observed by us. 
Correcting our results upward for the area covered by the 
Ha map in JOS increases the H2 luminosity given above by 
a factor of 10. With these corrections, we find that the H2 
luminosity obtained in this work is in good agreement with 
the JOS results. 

To conclude, upon correcting our results for the area 
covered by the total extent of the ionised gas nebulae, we 
find H2 luminosities i(H2, A2S97)~10*^ erg s"^ and i(H2, 
S1S9)~10''2 erg s"^ for our clusters. The analysis performed 
above relies heavily on the assumption that the H2 gas is in a 
single temperature phase. If there is colder H2 gas pres ent in 
these clusters, as expected from CO observations by lEdg j 
(|200ll ) and mid-infrared spectroscopy (Jaffe & Bremer in 
prep.) we underestimate the total H2 luminosity here. 

7.5 Mass of the Warm Molecular Gas 

In the previous section we found that a LTE model provides 
a good description of the observed H2 line ratios in A2S97 
and assumed t hat t he sa me is true for a S159. Fo ll owing 
IScoville et al] (| 19821 ) and IStorchi-Bergmann et all (|2009l i 
the mass of a H2 gas in LTE conditions with a single 
kinetic temperature, can be calculated using the following 
equation, 

2mpfH2i-osi47rD^ 
= 7 -. T— (3) 

= (S.08 X 10^^) ( ( (4) 
^ ' \^ergs-icm-2 J \ Mpc J ^ ' 
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Here m-p is the proton mass, h is the Planck constant 
and D is the distance to the cluster core. -Fnai-osi is 
the flux of the H2 1-0 S(l) line and /^^^ j^g is its pop- 
ulation fraction. j4h21-osi is the transition probability 
(^H2i-osi = 3.47x10"'^ s"^) and V is the frequency of the 
H2 1-0 S(l) line. The right-hand side of equation (2) is 
obtained by assuming a vibration temperature TVib = 2000 
K for the gas, in which case /^=i_j^3 = 1.22x10"^ In the 
above formula is given in units of solar masses. 

The resulting H2 gas masses for the regions investigated 
here are given in Tables iQland l 121 The total H2 gas mass inte- 
grated over all observed flelds is 4.5x10^ ± 4.2x10^ Mq for 
A2597 and 5.2x10^ ± 2.9x10^ Mq for S159. Since our ob- 
servations do not cover the entire cluster core these numbers 
underestimate the total amount of molecular gas present in 
these clusters. We furthermore note that there is evidence 
from mid-infrared spectroscopy from Jaffe & Bremer (in 
prep.) that a colder H2 phase ( T ~ 300 K) exists in A2597 
and S159. If this colder H2 gas coexists with the warmer 
H2 gas observed here this means that we strongly under- 
estimate t he H2 mass present in these clusters. CO obser- 
vations by lEdgd l|200ll ) indicate that even colder II2 gas is 
present in A2597 at a temperature of a few tens of kelvin. 
The total molecular gas mass inferred from CO for A2597 is 
8x10^ M0, although this value relies heavily on the CO to 
H2 conversion factor. To our knowledge there is no published 
CO detection for S159. 

7.6 Mass of the Ionised Gas 

Following IStorchi-Bergmann et al.l l|2009l ) the mass 
of the ionised hydrogen gas can be estimated as 
Mhii = mprieVmi- Here nip is the proton mass, no is 
the electron density and Vhii is the volume of the HII 
emitting region. Using j„a/ n^, = 9.788x10"^^ erg s "'^ crn *^ 
and = 0.339 (|Osterbrock fc Ferlandl l2006h . 

for rjc — 10^ cm~^ and T = 10* K), we can write 
Fp.c. = JPac ^Hiil>-' = 3.32xl0-2^x(n2ymi£|-'). Here 
Fpaa is the flux of the Pa a line, jp^^ is the volume 
emissivity coefficient of Pa a and D is the distance to the 
cluster core. The HII gas mass can be written as: 



Mhii = mpncVnii (5) 




Using the right-hand side of equation (3) the mass of 
the HII emitting gas, Mhii, is given in units of solar 
masses. iHeckman et all (| 19891 ) finds an electron density, 
Ue = 200 cm~^, for the central regions of A2597. Thus 
assuming an electron density Uc — 200 cm"'^, the mass of 
the HII emitting gas in the selected regions of Abell 2597 
and Sersic 159-03 can be calculated. 

The resulting HII gas masses are given in Tables [9] and 



1121 The total HII gas mass integrated over all observed 
fields is 9.7x10** ± 3.2x10^ M© for A2597 and 7.6x10^ ± 
8.4x10'' Mq for S159. Since our observations do not cover 
the entire cluster core these numbers underestimate the to- 
tal amount of ionised gas present in these clusters. The total 
ionised gas mass to warm molecular gas mass is similar in 
both clusters; Mhii/Mh2 ~ 10^ . The HII gas mass de rived 
here for A2597 agrees well with'Heckman et al.' ('1989') who 
find Mhii = 5.8x10® Mq for a region with radius 2.5 kpc 
centered on the BCG nucleus. 

7.7 Stability of the Filaments 

From our observations we see that the ionised and molec- 
ular gas is locked up in large scale filamentary structures 
surrounding the BCG. If the filaments in our the clusters 
are similar to t he ones observed i n more nearby clusters, 
such as Perseus l|Fabian et al.llioO^ '). then we do not resolve 
their true physical thickness. Furthermore if we assume that 
the filaments observed here are connected to the global dis- 
tribution of ionised gas in these clusters then its clear that 
these extent far beyond the distance observed here by us 
(J05; lO'Dea et al.ll2004l '). 

This interpretation is strengthened by higher resolu- 
tion images in fo r example Ha and Ly a emission (DOO; 
lO'Dea et al]|2004l ). These show the existence of narrow gas 
filaments coincidental with the larger scale structures ob- 
served here. Considering this, even if these large structures 
consist of many unresolved filaments it seems plausible that 
all of these have a similar origin. The currently favoured 
picture is that these structures are related to past and cur- 
rent AGN outfiows. Another explanation may be that these 
structures arise in regions of these clusters where heating 
and cooling rates are not balanced and thus that we observe 
residual cooling here. 

Not only is the origin of these structures unclear, there 
are many other puzzling aspects concerning them. On large 
scales these structures appear to connect all the way down 
to the nucleus of the BCG. One may speculate on whether 
this gas falls in towards the black hole thereby feeding it and 
sustaining its activity. However, the velocities observed for 
the fi laments studied here or in other clusters (|Hatch et al.l 
[2OO5I) does not agree with gas that is freely infalling from 
distances far away from the potential centre. We observe the 
gas to be moving at velocities of about 100 km s~^, whereas 
for freely infalling gas we would expect to observe velocities 
of about 600 km s"'^. 

We would expect tidal forces to rip these filaments 
apart as discussed in iFabian et al.l (|2008l 'l. These authors 
suggest that magnetic fields may help stabilise these fila- 
ments against these gravitational forces. For a typical fila- 
ment in Perseus, or one of our clusters, the required magnetic 
field is around a few hundred /iG From equipartition argu- 
ments the ICM typically has magnetic field strengths less 
than 10 jj,G. Larger magnetic field strengths, up to about 
50 /xG, have been found in the centres of clust ers but these 
are usuall y related t o current AGN outfiows ijTavlor et al.l 
[2007; Govoni fc Feret ti 2004) . Significantly higher magnetic 
field strengths in the filaments are necessary to stabilise 
them against gravitational forces. 

Some of the gaseous filaments appear, in projection, 
to be related to the current AGN outfiows, see Fig. [6] and 
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[9l These outflows maybe be responsible for pushing gas 
outwards, thereby in fact creating the observed filaments 
llHatch et "ai]l2005l ). This may explain why some fllaments do 
not show large infall velocities. Ifowever, there are also fila- 
ments that are not related to current radio outflows. If these 
fllaments were deposited at their current location by previ- 
ous AGN outflows, then this gas has had a long time to cool 
and should now be falling back at high velocities towards 
the potential centre. No signs of high velocity infalling gas 
is observed in either of our clusters. We will discuss the re- 
lation between the filaments and the radio emission in more 
detail below. 



8 X-RAY AND RADIO EMISSION 

In the above sections we have shown that there is a strong 
spatial and dynamical relation between the HII and H2 gas 
phases in the cores of A2597 and S159. Below we investigate 
how these gas phases relate to X-ray and Radio emission. 



8.1 X-ray emission 

We find that globally the H2 and HII line emission in A2597 
appears to match the X-ray emission as imaged by CHAN- 
DRA rather well, see Fig. [T] However, upon examining the 
detailed distribution. Fig. 1131 we find that the X-ray bright- 
ness profile is considerable shallower than that of the colder 
gas. 

In S159, see Fig. 1141 we observe that the peak X-ray 
emission is significantly offset from the BCG nucleus and the 
peak H2, HII emission. One may perhaps argue for a relation 
between the X-ray and Pa a emission along the northern 
filament, but Pa a is difficult to detect in the nuclear region. 
There appears to be no relation between Pa a and X-ray 
emission along the western filament in S159. This leads us 
to conclude that is no strong correlation between the X-ray 
and either the H2 or HII phase in A2597 and S159. 

The total X-ray luminosity is 

ix(A2597) = 2.1x10** erg s'^ (|de Grandil Il999l. 0.5- 
2.0 k eV) and Lx(S159) = 9.6x10*^ erg s"^ (|de Grandil 
I1999I . 0.5-2.0 keV) in A2597 and S159 respectively. Taking 
into account the difference in areas between this work and 
that of J05, the total warm H2 luminosity for these clusters 
is a factor 10-100 times less than their X-ray luminosities 
and not of the same order of magnitude as previously 
claimed in J05. Adding the colder H2 gas to this analysis 
(Jaffe & Bremer in prep.) to obtain a total warm plus cold 
H2 luminosity may change this conclusion. 



8.2 Radio emission 

In Figs. H [3 El [91 [To] and [m we show VLA 8.4 GHz radio 
contours on top of the H2 1-0 S(3) and Pa a surface bright- 
ness, velocity and velocity dispersion maps. The total radio 
power is 3.1x10*^ erg s~^ f or A2597 and 2.1x10*^ erg s"^ 
for S159 (|Birzan et akllioosl . integrating the radio spectrum 
between 10 MHz and 10 GHz). 

Our observations of the molecular and ionised gas in 
A2597 and S159 do not have the spatial resolution to inves- 
tigate the detailed sub-kpc scale correspondence between the 



AGN outflows and the gas. On kpc scales we note small en- 
hancements in the HII and H2 intensity along the lower parts 
of the projected radio lobes in A2597. Most notably along 
the northern lobe. DOO, using higher resolution HST imag- 
ing, have previously shown this in more detail for A2597. 
Similarly in S159 we observe that on kpc-scales some of the 
gas appears to lie along the radio lobes. In neither S159 and 
A2597 is there a clear correlation between the radio emission 
and the cold gas. 

It has been postulated that gas observed on the edges 
of radio lobes is gas th at was uplifted an d deposited here 
through AGN outflows ([Hatch et al.ll2005l 'l. Another possi- 
bility in which the AGN outflows may create the observed 
structures is via pressure driven compression of the in-situ 
thermal gas by the outward expanding non-thermal plasma. 
In this latter scenario the pressure driven compression in- 
creases the local gas density leading to an increase in the 
cooling rate here. 

That the gas observed along the radio lobes is a purely 
projectional effect and has nothing to do with the radio out- 
flow seems unlikely. Strong evidence that an AGN interacts 
with the gas in the nuclear region of a BCG comes, for exam- 
ple, from the flux enhancement and high velocity dispersion 
structure observed in A2597 here and in DOO. These struc- 
tures coincide with the current radio outflows. Many other 
observations, mostly done in X-rays, have shown that there 
is a global correlation between the radio outflows in BCGs 



and the hot X-ray gas sur r ounding it (e.g..[McNamara et al 
l200]j; [Fabian et all [2OO2I : [Clarke et all [2005[ : [Birzan et al 



20081). We do caution the reader that even though this global 
correlation is found in many cool core clusters that there 
are also cases where one observes similar amounts of gas 
in structures that appear to have no relation at all to the 
current radio outflows. One example of this is the western 
fllament in S159. 



8.2.1 Alignment of the radio and gas rotation axes ? 

In A2597 we note that, in projection, there is a rough align- 
ment between the radio axis and the axis of rotation for the 
HII and H2 gas. The same seems to be true for S159 although 
the evidence for gas rotation is less clear here. Three other 
clusters, AbeU 1664, PKS 0745- 19 and A2204, have r ecently 
been observed with SINFONI bv [Wilman et al.l ([20091 '). They 
show that the HII and H2 gas in PKS 0745- 1 9 rota tes. High 
resolution radio images bv [Baum fc O'Deal | 199l[ ) show us 
that again that the radio axis and gas rotation axis roughly 
agree. For the other two clusters the situation is unclear. 
The velocity fields of the gas in AbeU 1664 and Abell 2204 
do not show clear e vidence for gas rotation, and radio im- 
ages b y [Govoni et a l. (2001) and Sanders Fabian & Tavloil 
([2OO9I ) do not allow for a reliable identification of the jet 
axis. Higher spatial resolution spectroscopy and radio imag- 
ing will be required to investigate the possible alignment of 
the radio axis and the gas rotation axis in more detail. 



8.2.2 Kinetic energy of the outflows 

It has been proposed in the literature that AGN feedback 
in the form of radio outfiows, jets could deliver the required 
heat that keeps the gas in cluster cores from cooling (e.g.. 
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iBirzan et al] |2004| : iDunn fc FabianI |2006| ). The kinetic 
power in the molecular and ionised phases may be a useful 
indicator for estimating the total energy input into the 
cluster medium by the current jet. Here we will estimate 
the kinetic power from the kinetic energy of the disturbed 
gas and its turbulent dissipation. The kinetic energy Ek 
and kinetic (turbulent) power Pk may be written as; 



Pk 



2 

2r 



3Ma^ 
2 



(7) 
(8) 



Here v = VSa with the a the velocity dispersion of 
the gas and r is the thickness of the high dispersion 
features. We will perform this estimate using region Al 
in A2597. This region contains most of the gas whose 
dispersion is clearly related to the AGN as decribed in 
Section gl The mass of the HII gas in Al is 4.1x10*^ Mq 
and its velocity dispersion is 256 km s~^. We use only the 
HII gas here as this gas component dominates in mass over 
the warm H2 gas. We thus find that a total kinetic energy 
of Sk(A2597:A1) = 8.1x10^'* erg and a kinetic energy 
power Pk (A2597:A1) = 5.8x10*° erg s"\ where we used 
r — 2 kpc. 

For S159 this calculation is more uncertain due to the 
lower significance of the Pa a line in the nuclear region Bl. 
If we use that the mass of the HII gas in Bl is 5.8x10* 
M0 and its velocity dispersion is 127 km s~^. We then find 
that a total kinetic energy of £;k(S159:B1) = 2.8x10^^ erg 
and a kinetic power Pk(S159:B1) = 1.0x10^* erg s"^ is in- 
jected into stirring up the gas, again using r — 2 kpc. These 
numbers are rather low for S159. If we use the H2 1-0 S(3) 
velocity dispersion for the Pa a line and increase the Pa a 
flux such that the Pa a/H2 1-0 S(3) ratio is the same as in 
region Al for A2597 we find that Ek and Pk both increase 
by about a factor of 10 for S159. 

We compare this mechanical power injected into the 
cool gas to the current radio power. We assume that the 
stirring of the gas in the nuclei of A2597 and S159 is related 
to the most recent AGN outburst and thus we look only at 
the recent radio emission. This is done to avoid the older 
radio plasma, from previous outbursts, which will dominate 
the radio power due to its emissivity at low frequencies. 

The 8.4 GHz radio emission is a good indicator of re- 
cent AGN activity. The synchrotron lifetime of electrons 
emitting at 8.4 GHz is about 10^~^ years, for a mag- 
netic field strength of 30 fiG (|Tavlor et all |2007| 1. The 
8.4 GHz specific luminosities Li, in regions Al and Bl 
are 1.75x10^^ erg s'^ Hz'^ and 3.88x10^^* erg s~^ Hz"\ 
One can then estimate the total current radio power Prad 
in these regions by Prad — u x L^, with u = 8.4 GHz. 
We find Prad(A2597:Al) = 1.5x10*^ erg s"^ and 
Prad(S159:Bl) = 3.3x10^^ erg s'V 

The kinetic power inserted into the cool gas is thus 
smaller than the radio power in A2597 by a factor 2-3. For 
S159 we find that Pk is smaller than Prad by more than an 
order of magnitude. If we allow for the higher dispersion and 
flux as mentioned above then this ratio becomes similar to 



that observed in A2597. Adding the warm ( T ~ 2300 K) H2 
gas to the above analysis does not significantly change the 
derived kinetic powers due to its low mass. We thus find that 
within the nuclear regions the current radio jets in A2597 
and S159 are radiating efficiently with respect to the warm 
HII and H2 gas. 

The kinetic power deduced here should though be seen 
as a lower limit to the total kinetic power of the cool 
(T^^IO* K) gas in these clusters since the colder gas phases 
have not been taken into account. 094 find a kinetic energy 
Ek ~10^® erg and a kinetic power Pk ~10*^ ergs~^ for the 
HI gas in a region with radius 2.5 kpc centered on the BCG 
nucleus in A2597. Similarly the cold (T^IO^ K) molecular 
gas may contain a significant amount of kinetic energy. Stud- 
ies of X-ray cavities find evidence for inefficiently radiating 
jets in that the mechanical power of the jet often can be 
orders of magnitude larger than the radio power. However, 
when looking at X-ray cavities one may be looking back in 
time towards an epoch in which the balance between the me- 
chanical and radiative power of the radio source is different 
from the current one. 



8.2.3 Low velocity shocks 

The argument that AGN feedback is responsible for reheat- 
ing the ICM gas in cool core clusters is made purely on the 
global energetics of the problem and does not consider the 
microphysics of the energy transfer between the non-thermal 
radio plasma and thermal ICM plasma, nor does it con- 
sider the various scales involved in this energy transfer. Fast 
shocks {v J5 40 km s~^) can be ruled out based on the ab- 
sence of typical shock tracers like the [O HI] 4363 Angstrom 
line in optical spectra (VD97) and the high H2/HII ratios 
(|Donahue et al.ll2000l ). However, another way of delivering 
the required heat for the cooler gas observed here is via 
low velocity shocks in a dense gas. Such a model was re- 
cently invoked to exp lain the strong H2 em ission observed in 
'Stephan's Qui ntet' (iGuiUard et aklbOOgl'). From X-ray ob - 
servations (e.g. 'McNa mara et aiT l200l': 'Fa bian et al.ll2006l ) 
and jet models (iSutherland fc Bicknell 20071 e.g.) it is found 
that AGN outflows may give rise to low velocity C-type 
shocks through a turbulent cascade of the expanding jet mo- 
mentum in a inhomogeneous gas. 

We have investigat ed lo w velocity C-type shocks by 
iKristensen et al.l |20o3) and iFlower fc Pineau des Forets I 
(|2003t ) 1,0 model the molecular ratios in region Al for 
A2597. The following parameter range was allowed for, 
shock velocity vs = [15-30] km s~^, pre-shock density 
riH = [10^-10^] cm~^ and magnetic field strength B = [0.05- 
30] mG. The best-fitting model to the molecular line ratios 
is vs = 22 km s~\ nn = 5x10*^ cm"^ and B = 5.6 mG. 
This model is shown by the purple line in Fig. [TU] It implies 
a pre-shock density above the critical density and thus that 
the gas is already in LTE prior to the shock. This model 
furthermore implies a very high magnetic field in the molec- 
ular gas, something which has not yet been observed in our 
objects. 
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9 SUMMARY 

Warm molecular gas is observed in the central 20 kpc of 
the A2597 and S159 galaxy clusters. All of the H2 1-0 S- 
transitions and H2 2-1 S-transitions redshifted into the K- 
band and none of the higher ro-vibrational H2 transitions 
are detected in A2597. Ionised gas is observed in A2597 
through the detection of the Pa a, Br 7, Br 5 and Fe II 
(1.81 /xm) lines. For S159 only two H2 1-0 S-transitions are 
detected. There is tentative evidence for the presence of the 
other II2 1-0 S-transitions, but these are strongly affected 
by atmospheric effects. None of the higher ro-vibrational H2 
lines were detected. Ionised gas is observed in S159 through 
the detection of the Pa a line. Higher ionisation lines, such 
as [Si IV], are not detected in either S159 or A2597. 

Molecular and ionised gas is detected out to 22.5 kpc 
from the nucleus of PGC071390 in A2597 and 18.0 kpc from 
the nucleus of ESO291-G009 in S159. The fields observed 
by us do not cover the entire extent over which molecular 
and ionised emission has previously been detected in these 
clusters. Our observations thus miss a significant fraction of 
the total emission and gas mass. Based on the Ha maps in 
JOS we estimate that we have observed about one tenth of 
the total HII and H2 gas present in A2597 and S159. 

The gas morphology in A2597 is more diffuse than in 
S159. In both clusters the ionised and molecular gas is ob- 
served to be distributed in clumps and filamentary struc- 
tures surrounding their BCGs. This observation brings these 
higher redshift cool core clusters closer to the situation ob- 
served in more nearby clusters, such as Perseus and Centau- 
rus, that show a wealth of small scale filamentary structure 
(e.g. Fabian ct al. 2008; Crawf ord et al. 200 5). 

In both A2597 and S159 the ionised and molecular gas 
phases trace each other closely on kpc-scales, in both ex- 
tent and dynamical state. The H2 line ratios and the H2 
to HII line ratios only vary slowly as a function of posi- 
tion throughout the investigated regions, with the excep- 
tion of the nuclear regions. This implies that there is a tight 
coupling between these gas phases and that there may be 
a single excitation mechanism responsible for the observed 
emission. 

In both clusters the Pa a to H2 ratio is found to be 
lower within their nuclear regions than outside. In A2597 
this decrease is small as the Pa q:/H2 1-0 S(3) ratio decreases 
from 0.94±0.04 in region A2 to 0.77±0.07 in the nuclear 
region Al. In S159 the decrease is much larger and we find 
that the Pa a/H2 1-0 S(3) ratio decreases from 0.66±0.12 in 
region B3 to 0.19±0.08 in the nuclear region Bl. This may 
imply that the physical conditions differ within the nuclear 
regions as compared to further out. 

The velocity structure of the ionised and molecular gas 
within the nuclear regions of A2597 and S159 is consistent 
with gas rotating around the nucleus. The data suggests a 
possible alignment between the axis of rotation and the radio 
axis. Further high spatial resolution, spectroscopic observa- 
tions are required to disentangle this rotation from velocity 
gradients observed along the filamentary gas structures and 
to confirm the axis of rotation. 

The high velocity dispersion of the gas within the nu- 
clear regions clearly shows that the AGNs harboured by 
these BCGs interact with the molecular and ionised gas sit- 
uated here. This is especially obvious in A2597 where the 



observed increase in the velocity dispersion is suggestive of 
a turbulent wake marking the passage of the radio lobes 
through the surrouding medium, or of direct acceleration of 
the dense material by the AGN jet. 

We estimate that the mechanical power of the AGN flow 
is somewhat smaller than its radiative radio power within 
the nuclear regions of A2597 and S159 with respect to the 
ionised and warm molecular gas. If we were to include the 
mechanical energy contained within the colder atomic and 
molecular gas phases than this situation will reverse. 

The velocities observed for the filamentary structures 
are less than 150 km s~^. This does not agree with freely 
infalling gas which should have much higher velocities, i.e. 
~600 km s~^. The filaments do show smooth velocity gradi- 
ents. Possible explanations for this are; (i) we see a single fil- 
ament that is being stretched and bended by its proper mo- 
tion and gravitational forces, (ii) we see a blend of a number 
of unresolved filaments along the line of sight. From observa- 
tions of nearby clusters suc h as Perseus the latter explana- 
tion a ppears most likely fe.g. lFabian et a02008l : lHatch et al.l 
l2005h . 

The velocity dispersion of the ionised and molecular 
gas decreases with distance to the nucleus. It drops below 
100 km s~^ at distances of a few kpc from the nucleus which 
implies that the gas is kinematically cold outside the nu- 
clear region. The decrease in velocity dispersion is strongest 
within the first ~2 kpc from the nucleus and slows down af- 
ter that. This point may mark a transition in where the AGN 
is able to strongly interact with its surrounding medium. 

The low velocity and velocity dispersion of the gas in 
the filamentary structures indicates that the gas here resides 
in well defined structures. The gas may still be clumpy but 
individual clumps must have a spread in velocity no larger 
than the derived dispersion. It also implies that these struc- 
tures need some form of support, which can not be kinemat- 
ical, to keep them from plummeting towards the potential 
centre. A pos sible source o f support are magnetic fields as 
suggested bv I Fabian et al.l ((2OO8). This requires magnetic 
field strengths of a few hundred /iG. So far there is no ob- 
servational evidence for the existence of these magnetic fields 
in galaxy clusters. 

The molecular and ionised gas emission observed in 
A2597 and S159 does not match the observed X-ray emis- 
sion very well. In S159 there is no match in the observed 
brightness profiles. In A2597 there is a rough global agree- 
ment in the peak emission but the detailed X-ray brightness 
profiles are significantly shallower than observed for the H2 
and HII gas. 

There is no evidence for a strong correlation between 
either H2 or HII and the radio emission in A2597 and S159. 
If the AGN is main source of heating and excitation for 
these gas phases, then we would expect to find this gas near 
the current AGN outflows at all times. This is especially 
important for the H2 gas since it has a cooling time of order 
years. One prime example where we find H2 gas that appears 
to have no relation to the current outflows is the western 
filament in S159. The emission line spectra for all regions in 
S159 are similar. This then implies that the physical proces 
that is heating the gas in these regions is similar and likely 
not directly related to the current AGN outflows. 

In A2597 the H2 line emitting gas is well described by 
an LTE model with a single excitation temperature Texc. 
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We find that Tcxc is about 2300 K for tiie different regions 
investigated. We find sligiitiy iiiglier excitation temperatures 
for the molecular gas in the filaments as compared to the 
nuclear region. 

In S159 we have only been able to reliably measure two 
H2 lines. The temperature and excitation of the H2 gas in 
this cluster thus remains largely unconstrained. Assuming a 
single temperature LTE model, we find that Texc is about 
2300 K for the different regions investigated. The temper- 
ature derived in this way for the nuclear region in S159 is 
above the dissociation temperature of H2 and considered 
unphysical. No clear trends in temperature with distance to 
the nucleus are observed for S159. 

If the molecular gas in A2597 and S159 is in LTE then 
this implies that the H2 gas is not in pressure equilibrium 
with the HII ff £LS , £LS has previously been pointed out by JOS. 
We find that Mhii/Mh2 ~ 10^ in A2597 and S159. We note 
that Mhj is likely to be seriously underestimated since only 
the K-band H2 lines tracing the warm molecular gas have 
been considered here. 

For a region with radius 2.5 kpc centered on the nucleus 
of the A2597 BCG we can make an inventory of the cool gas. 
In agreement with lHeckman et al.l (| 19891 ) we find an ionised 
gas mass Mhii -6x10® Mq. We find a warm molecular gas 
mass MH^.warm -3x10* Mq. lO'Dea et all l|l994) find an 
neutral hydrogen mass Mhi — 10^ M©. Using much lower 
spatial resolution Edg3 l|200ll 'l find a spatially unresolved 
cold molecular gas mass MH2,coid=8x 10^ Mq in the central 
22" of A2597. 

We are able to fit the II2 line emitting gas by a low 
velocity shock model. Such a model requires high pre-shock 
densities and very high magnetic fields. Low velocity shocks 
may be a viable option for heating the H2 gas in cool core 
clusters, but this requires further modelling and observa- 
tions. 

Using recent Spitzer observations we find that another 
lower temperature II2 gas component exists at about 300 K 
in A2597 and S159 (Jaffe & Bremer in prep.). This shows 
that a single temperature LTE model is not a good descrip- 
tion for all of the H2 gas present in cool core clusters and 
again brings the situation in our clusters c loser to that ob- 
serve d in more nearby cool core clusters l|johnstone et al.] 
I2OO7I ). This is also implies that we significantly underesti- 
mate the total H2 luminosity and gas mass by using only 
the H2 K-band spectrum. 

Our K-band spectra of cool core clusters do not re- 
semble typical K-band spectra of other line emitting, ex- 
tragalactic objects such as Ultraluminous Infrared Galaxies 
(ULIRG) and AGNs that are not in BCGs. For comparisons 
we can look at the Br 7 and II2 1-0 S(l) line. In ULIRGs 
theses line have a similar intensity and width and in AGNs 
that are not in BCGs the Br & line usually has a higher inten- 
sity and lar ger width than the H2 1-0 S(l) line (IDavies et al.l 
l2003l . 120051 ). In our spectra for AbeU 2597 we find that these 
lines have similar widths, but the II2 1-0 S(l) line is more 
than ten times brighter than the Br 7 line. 



10 CONCLUSIONS 

Above we have presented the first K-band integral-field spec- 
troscopic observations of the extended molecular and ionised 



gas distributions in the core regions of the A2597 and S159 
galax;y clusters. These observations add one extra dimension 
to the analysis of the molecular and ionised gas in these clus- 
ters as compared to our previous investigation in J05. This 
allows us to study the distribution, kinematics and thermal 
structure of this gas over a much larger area than in our 
previous investigation. The spatial resolution of our obser- 
vations limit us to studying this gas on kpc-scales. Below we 
summarise our conclusions. 



• Line emission from molecular and ionised gas is found 
in all observed fields. We confirm the conclusion in J05 that 
molecular and ionised gas are tightly coupled out to 20 kpc 
from the BCG nucleus in A2597 and S159. 

• The molecular and ionised gas is distributed in filamen- 
tary structures surrounding the BCG. The gas morphology 
is more diffuse in A2597 than in S159. 

• In all regions where we have sufficient signal to noise we 
find that the H2 gas can be described an by an LTE model 
with a temperature T ^ 2300 K. LTE implies a high-density 
for this gas. 

• Kinematically there is a fairly clear separation between 
the central few kpc, energetically dominated by the AGN, 
and an area outside where the support of the HII and II2 
gas remains to be explained. 

• Within the nuclear region of A2597 and S159 the kinetic 
luminosity of the HII and warm H2 gas is somewhat smaller 
than the current radio luminosity. 

• The very high dispersion regions observed in H2 and 
HII may represent a turbulent wake of the relative motion 
of galaxy and surrounding medium, or the interaction of the 
current jet with this medium. 

• A2597 shows a better correlation between the morphol- 
ogy of the cool gas and X-ray emission than S159. The de- 
tailed radial brightness profiles of the X-ray and either the 
H2 or HII gas do not match in either cluster. 

• The high frequency radio emission and cool gas have 
the same overall scale size. There is no strong correlation 
between either the H2 or HII gas and the radio structures, 
but there are weak enhancements in the surface brightness 
of this gas along the northern radio lobe in A2597. 

• The data suggests that the HII and H2 gas situated 
within a few kpc from the nucleus rotates around an axis 
parallel to the radio jet, this needs further confirmation. 



A more detailed investigation into the heating and cooling 
of the intracluster medium will need to take into account the 
microphysics of the interaction between the various phases 
situated within the intracluster medium; the hot X-ray gas, 
the cool HII, H2 gas and the radio emitting plasma. Once 
such a model is in place, the next step will be to compare 
detailed maps of local heating versus local cooling in cool 
core clusters. We will follow this work up in the near future 
with optical integral field spectroscopy (Oonk et al. in prep.) 
and Spitzer IRS spectroscopy (Jaffe & Bremer in prep.). A 
more detailed investigation of the excitation mechanisms for 
the molecular hydrogen gas will also be presented in a future 
paper. 
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Cluster 


2 


Dl [Mpc] 


[kpc arcsec ^] 


Lx [10** erg s-i] 


ii.4GHz [10=*! erg s-i Hz-i] 


t-ohs [h] 


Sersic 159-03 


0.0564 


245 


1.1 


2.1 


29.3 


4.2 


Abell 2597 


0.0821 


363 


1.5 


1.0 


1.6 


8.5 



Table 1. Targets. Columns 1, 2 give the cluster name and the redshift of its BCG. In column 3 we give the luminosity distance to the 
BCG in units of Mpc. In column 4 we give t he angular size sc ale at the distance of the BCG. In column 5 we give the X-ray luminosity 
of the cluster in th e 0.5-2.0 keV band from lde Grandil lll999l ') and in column 6 we give the radio power of the cluster at 1.4 Ghz from 
iBirzan et al] ll2008h In column 7 we give the on-source integration time in units of hours. An equal amount of time was spend off-source 
observing the sky. 



Resolution 


no smoothing 


ss20 


ss22 


ss44 


Spatial [arcsec] 


0.9 


1.0 


1.0 


1.1 


Spectral [km s~^] 


38 


38 


51 


73 



Table 2. Spatial and spectral Resolution. Dependence of the spatial and spectral resolution on smoothing. The spatial resolution is given 
in units of arcsec and the spectral resolution in terms of km s~^. Column 1 gives the type of resolution. Column 2 gives the resolution for 
no smoothing. Column 3 gives the resolution after smoothing with a two pixel FWHM Gaussian in the spatial domain. Column 4 gives 
the resolution after smoothing with two pixel FWHM Gaussians in both the spatial and spectral domain. Column 5 gives the resolution 
after smoothing with four pixel FWHM Gaussians in both the spatial and spectra domain. 



Field [10 erg s ^ cm ^ /^m ^] 


no smoothing 


ss20 


ss44 


Central 


24.9 


12.6 


3.7 


South-East 


40.0 


19.8 


5.4 


South-West 


25.8 


12.9 


3.5 


North 


30.1 


15.0 


4.0 



Table 3. ABELL 2597 Sensitivity. The RMS noise is given in units of 10~^^ erg s~^ cm~^ /im~^ for an area of 0.125x0.125 arcsec^ 
which is equivalent to one spatial pixel. The RMS noise was calculated in the interval 2.07-2.30 /^m with a sampling of 2.45x10"'* ^m 
which is equivalent to one spectral pixel. Column 1 gives the observed field. Column 2 gives the RMS noise for no smoothing. Column 
3 gives the RMS noise after smoothing with a two pixel FWHM Gaussian in the spatial domain. Column 4 gives the RMS noise after 
smoothing with a four pixel FWHM Gaussian in both the spatial and spectral domain. The RMS noise does not vary much across the 
observed field. 



Field [10 erg s ^ cm ^ /^m *] 


no smoothing 


ss22 


ss44 


South-East 


33.1 


13.2 


9.0 


South-West 


30.5 


10.2 


4.3 


North 


36.7 


12.2 


5.0 



Table 4. SERSIC 159-03 Sensitivity. The RMS noise is given in units of 10"^'^ erg 8"^ cm'^ /im"! for an area of 0.125x0.125 arcsec^ 
which is equivalent to one spatial pixel. The RMS noise was calculated in the interval 2.07-2.30 /xm with a sampling of 2.45 xlO"'* /^m 
which is equivalent to one spectral pixel. Column 1 gives the observed field. Column 2 gives the RMS noise for no smoothing. Column 
3 gives the RMS noise after smoothing with two pixel FWHM Gaussians in both the spatial and spectral domain. Column 4 gives the 
RMS noise after smoothing with four pixel FWHM Gaussian in both the spatial and spectral domain. The RMS noise does not vary 
much across the observed field. 
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Line Name 


Flux [10 erg s ^ cm ^] 


Luminosity [10 


H2 1-0 b(0) (2.2235 ^va)^ 


121.7±11.3 


19.2±1.8 


±12 i-U o(-Lj (Z.izio fiTH) 


0Dy.Uiti4.D 


lUo.oitz.o 


H2 1-0 S(2) (2.0338 Aim) 


280.2±8.4 


44.2±1.3 


H2 1-0 S(3) (1.9576 Aim) 


732.0±19.3 


115.4±3.0 


H2 1-0 S(4) (1.8920 Aim) 


129.7±7.5 


20.5±1.2 


H2 1-0 S(5) (1.8358 Aim) 


462.2±11.7 


72.9±1.8 


H2 2-1 S(2) (2.1542 Aim)2 


22.0±4.1 


3.5±0.6 


H2 2-1 S(3) (2.0735 Aim)i 


38.0±2.2 


6.0±0.3 


H2 2-1 S(5) (1.9449 Aim),Br <5 (1.9451 fim)'' 


39.3±4.2 


6.2±0.7 


Br 7 (2.1661 Atm)2 


30.1±4.6 


4.7±0.7 


Pa a (1.8756 Atm) 


613.0±20.3 


96.7±3.2 


Fe II (1.8100 Atm)l 


86.5±5.6 


13.6±0.9 



Table 5. ABELL 2597 Integrated line fluxes and luminosities. The integrated line fluxes are obtained by collapsing the data cube into a 
single spectrum and fitting the lines by a single Gaussian. The data was smoothed by a four pixel FWHM in both the spatial and spectral 
domain. The results were inspected by eye and the errors were estimated using Monte-Carlo simulations. The spectrum also hints at the 
presence of the H2 2-1 S(4) line, but a reliable flux could not be derived for this line. Fluxes are given in units of 10~^^ erg s~^ cm~^ 
and luminosities are given in units of 10^^ erg s~^. 

1 The H2 1-0 S(0), H2 2-1 S(3), H2 2-1 S(5), Br & and Fe II Unes have been integrated over the central field only. 
^ The H2 2-1 S(2) and Br 7 lines have been integrated over the combined areas Al and A2 only, see Fig.|3] 



Line Name 


Flux [10 ^"^ erg s cm 


Luminosity [10^® erg s ^] 


H2 1-0 S(l) (2.1218 Atm) 


171.0±11.9 


12.3±0.9 


H2 1-0 S(3) (1.9576 Aim) 


217.9±18.1 


15.7±1.3 


Pa a (1.8756 Aim) 


104.4±11.6 


7.5±0.8 



Table 6. SERSIC 159-03 Integrated line fluxes and luminosities. The integrated line fluxes are obtained by collapsing the data cube 
into a single spectrum and fltting the lines by a single Gaussian. The data was smoothed by a four pixel FWHM in both the spatial and 
spectral domain. The results were inspected by eye and the errors were estimated using Monte-Carlo simulations. The spectrum also 
hints at the presence of the H2 1-0 S(0), 1-0 S(2) and 1-0 S(4) lines, but reliable fluxes could not be derived for these lines. Fluxes are 
given in units of 10~^^ erg s~^ cm"'^ and luminosities are given in units of 10^^ erg s"'^. 



Line/Area [10-1'' erg s-i cm-2] Al A2 A3 A4 A5 A6 A7 



H2 1-0 S(0) 


69.0±8.5 


18.6±3.0 














H2 1-0 S(l) 


332.0±6.6 


84.1±3.0 


4.9±1 


2 


21.8±6.1 


4.6±1.1 


12.0±0.9 


3.2±1.1 


H2 1-0 S(2) 


129.0±4.1 


31.0±2.1 


(2.1) 




14.6±1.7 








H2 1-0 S(3) 


335.3±9.2 


85.2±3.0 


6. Oil 





32.6±4.3 


4.8±0.7 


10.2±1.3 


2.6±0.9 


H2 1-0 S(4) 


69.7±4.0 


21.5±2.0 


1.9±0 


9 


10.9±2.1 








H2 1-0 S(5) 


223.1±5.4 


53.0±2.5 


3.7±1 


1 


25.4±4.5 


3.5±0.7 




1.5±0.7 


H2 2-1 S(2) 


16.2±2.4 


5.5±2.2 














H2 2-1 S(3) 


26.9±3.0 


7.5±1.0 














H2 2-1 S(5) 


13.7±4.1 


2.8±1.0 














Br <5 


14.1±4.1 


4.3±1.0 














Br 7 


21.1±2.5 


6.5±2.1 














Pa a 


257.9±8.8 


79.7±3.6 


5.3±1 


2 


29.0±7.0 


3.6±0.7 




2.4±0.8 


Fe II 


50.2±7.2 


5.0±2.1 















Table 7. ABELL 2597 Integrated line fluxes for regions A1-A7. The integrated line fluxes are obtained by collapsing the data cube, 
within the selected region, into a single spectrum and fltting the lines by a single Gaussian. The smoothing performed for a given line is 
specified in the caption of the line spectra shown in Appendix [C] The results were inspected by eye and the errors were estimated using 
Monte-Carlo simulations. If Br 7 is detected the H2 2-1 S(5), Br 5 complex is disentangled by assuming that the Br 7/Br <5 ratio is 1.5 
(Osterbrock k. Ferland 2006). The uncertain H2 1-0 S(2) flux in region A3 could not be fltted well. The H2 1-0 S(2) flux given for A3 is 
based on flxing the width of this line to that of the H2 1-0 S(4) line. Fluxes are given in units of 10" "'^'^ erg s" cm" ^. 
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Line/ Area [km s ^] 


Al 


A2 


A3 


A4 


A5 


A6 


A7 


v{H2 1-0 S(3)) 


+2.9±6.1 


-29.3±5.7 


-f36.0±19.1 


-10.2±12.9 


-fl9.8±14.2 


-fl74.0±16.7 


-l-156.1±24.5 


(7(H2 1-0 S(3)) 


214.3±6.0 


145.8±6.1 


97.6±19.9 


88.8±13.4 


103.4±14.9 


119.5±17.4 


60.1±14.9 


v(Pa a) 


-21.2±7.5 


-59.2±8.3 


-(-58.9±26.7 


-3.8±29.4 


+21.6±21.8 




-l-152.4±51.8 


(T(Pa a) 


255.9±8.7 


176.9±8.3 


108.1±27.5 


118.2±29.1 


115.9±23.1 




109.0±56.0 



Table 8. ABELL 2597 Kinematics for regions A1-A7. For each line the top row gives the velocity, with respect to the systemic velocity, 
and the bottom row gives the velocity dispersion. These are derived from a single Gaussian line fit to the collapsed spectrum for a selected 
region. Here we present results only for the H2 1-0 S(3) and Pa a lines. All other lines follow the behaviour observed in these two lines 
within errors. The smoothing performed for a given line is specified in the caption of the line spectra shown in Appendix [C] The results 
were inspected by eye and the errors were estimated using Monte-Carlo simulations. Velocity and velocity dispersion are both given in 
units of km . 



Lino/ Area 


Al 


A2 


A3 


A4 


A5 


A6 


A7 


Texc,H2 [K] 

Mhj [102 M0] 
Mhii [104 Mq] 


2229±38 

222.2±4.4 

410.0±14.0 


2141±59 
56.3±2.0 
126.7±5.7 


2792±808 

3.3±0.8 

8.4±1.9 


2660±481 

14.6±4.1 

46.1±11.1 


2678±640 

3.1±0.7 

5.7±1.1 


1586±264 
8.0±0.6 


1863±592 

2.1±0.6 

3.8±1.3 



Table 9. ABELL 2597 Gas temperatures and masses for regions A1-A7. For each region an excitation temperature, Tcxc,H2 in units of 
Kelvin, is calculated for the H2 gas. Molecular and ionised gas masses are calculated using the equations in the text. The molecular gas 
mass Mhj is given in units of lO'^ Mq. The ionised gas mass Mhii is given in units of 10* Mq, using Ue = 200 cm"'' and Thii = 10* K. 
The temperature and mass calculations for the H2 gas assume LTE conditions. 



Line/Area [10 erg s ^ cm 


-2] Bl 


B2 


B3 


B4 


B5 


B6 


B7 


H2 1-0 S(l) 
H2 1-0 S(3) 
Pa a 


26.6±2.6 
42.2±5.3 
8.0±2.4 


20.4±1.8 
24.1±2.0 
12.9±1.5 


16.4±1.7 
18.8±1.4 
12.4±1.3 


3.1±0.9 
3.5±0.7 
3.1±0.8 


35.7±2.6 
40.2±2.8 
25.2±3.2 


12.2±5.6 
12.6±3.3 
11.4±5.4 


6.0±1.2 
7.1±0.9 
5.6±1.2 



Table 10. SERSIC 159-03 Integrated line fluxes for regions B1-B7. The integrated line fluxes are obtained by collapsing the data cube, 
within the selected region, into a single spectrum and fitting the lines by a single Gaussian. The smoothing performed for a given line is 
specified in the caption of the line spectra shown in Appendix [C] The results were inspected by eye and the errors were estimated using 
Monte-Carlo simulations. Fluxes are given in units of 10~^^ erg s~^ cm"^. 



Line/ Area [km s 


Bl 


B2 


B3 


B4 


B5 


B6 


B7 


v(H2 1-0 S(3)) 


-26.3±28.6 


+109.1±12.2 


+101.8±8.3 


-48.7±17.1 


-23.4±7.8 


-62.9±20.5 


-19.1±9.0 


(t(H2 1-0 S{3)) 


249.9±31.1 


132.2±12.6 


103.1±8.3 


77.8±17.6 


100.0±7.8 


65.1±14.3 


61.6±4.4 


v(Pa a) 


-6.7±39.6 


+120.9±15.0 


+105.6±9.6 


-41.8±34.7 


-27.6±13.1 


-58.1±48.8 


-19.8±20.6 


(T{Pa a) 


126.5±39.8 


120.3±15.5 


95.1±9.6 


121.6±36.1 


103.5±15.1 


97.9±53.3 


85.7±20.3 



Table 11. SERSIC 159-03 Kinematics for regions B1-B7. For each line the top row gives the velocity, with respect to the systemic 
velocity, and the bottom row gives the velocity dispersion. These are derived from a single Gaussian line fit to the collapsed spectrum 
for a selected region. Here we present results only for the H2 1-0 S(3) and Pa a lines. All other lines follow the behaviour observed in 
these two lines within errors. The smoothing performed for a given line is specified in the caption of the line spectra shown in Appendix 
[C] The results were inspected by eye and the errors were estimated using Monte-Carlo simulations. Velocity and velocity dispersion are 
both given in units of km . 



Line/ Area 


Bl 


B2 


B3 


B4 


B5 


B6 


B7 


Texc,H2 [K] 
Mh2 [102 Mq] 
Mhii [10* Mq] 


5655±3632 

8.0±0.8 

5.8±1.7 


2520±546 

6.2±0.5 

9.3±1.1 


2391±518 

5.0±0.5 

9.0±0.9 


2331±1359 

0.9±0.3 

2.2±0.6 


2321±385 

10.9±0.8 

18.2±2.3 


2032±1550 

3.7±1.7 

8.3±3.9 


2527±1073 

1.8±0.4 

4.1±0.9 



Table 12. SERSIC 159-03 H2 Gas temperatures and masses for regions B1-B7. For each region an excitation temperature, Toxc,H2 
in units of Kelvin, is calculated for the H2 gas. Molecular and ionised gas masses are calculated using the equations in the text. The 
molecular gas mass Mh2 is given in units of 10^ Mq. The ionised gas mass Mhii is given in units of 10* Mq, using rie = 200 cm~^ and 
^HII = 10* K. The temperature and mass calculations for the H2 gas assume LTE conditions. 
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Figure 1. ABELL 2597 Observed fields, {left) extent of the X-ray emission in the cluster core as observed by CHANDRA (exposure 
time 153.7 ks). VLT FORS ij-band contours are overlayed in black (solid line), {right) zoom in of the region marked by the black dashed 
box in the left image. The fields observed by SINFONI are indicated with a dashed line in the right image. H2 1-0 S(3) contours are 
overlayed in black (solid line), drawn starting at 1.5 X lO"^'^ erg s~^ cm~^ arcsec"^ in steps of 2° with n=0,l,2,... The blue plus sign 
marks the position of the stellar nucleus. 




Figure 2. SERSIC 159-03 Observed fields . {left) the extent of the X-ray emission in the cluster core as observed by CHANDRA 
(exposure time 10.1 ks), has been smoothed slightly. VLT FORS i?-band contours are overlayed in black (solid line), {right) zoom in of 
the region marked by the black dashed box in the left image. The fields observed by SINFONI are indicated with a dashed line in the 
right image. H2 1-0 S(3) contours are overlayed in black (solid line), drawn starting at 1.5 X 10~^^ erg s~^ cm"'^ arcsec"^ in steps of 
2° with n=0,l,2,... The blue plus sign marks the position of the stellar nucleus. 
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Figure 3. Selected regions in ABELL 2597 and SERSIC 159-03. The regions are sliown overlayed on SINFONI H2 1-0 S(3) surface 
brightness maps. The left image shows A2597 and the right image S159. The fields observed by SINFONI are indicated with a dashed 
black line. The dotted black lines show the regions used for obtaining the off-source spatial median in these fields. The solid black lines 
show the selected regions A1-A7 for A2597 and B1-B7 for S159, see Sections [73] and [721 The small green and red squares show the 
regions investigated along the filaments, see Sections 14.31 and 16.31 The blue plus sign marks the position of the stellar nucleus. 
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Figure 4. ABELL 2597 full K-band spectra. As an example of the quality of our SINFONI spectra we show full K-band spectra for 
two regions. The top image shows the spectrum for region Al and the bottom image shows the spectrum for region A4, see Fig. |3] and 
Section l7.ll These spectra were obtained after smoothing the data by four pixels in both the spatial and spectral planes. The horizontal 
axis is given in units of /im and the vertical axis is given in units of lO"'^* erg cm~^ s"'^ fim~^ . 
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Figure 5. SERSIC159-03 full K-band spectra. As an example of the quality of our SINFONI spectra we show full K-band spectra for 
two regions. The top image shows the spectrum for the combined regions Bl and B2. The bottom image shows the spectrum for region 
B5, see Fig. [3] and Section 17.21 These spectra were obtained after smoothing the data by four pixels in both the spatial and spectral 
planes. The horizontal axis is given in units of fivci and the vertical axis is given in units of 10"^^ erg cm^^ s^^ /xm^^. 
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Figure 6. ABELL 2597 Surface brightness. We stiow surface brightness maps for the H2 1-0 S(3) (top) and Pa a (bottom) lines in units 
of 10~^^ erg cm^^ arcsec^^. The maps are obtained by fitting the spectrum for each spatial pixel in the data cube by a single 
Gaussian. The data was smoothed by four pixels in both the s patial and spectral planes. The stellar nucleus is indicated by the black 
cross. VLA 8.4 GHz Radio Continuum contours (solid black line. lSarazin et aLlll995h and CHANDRA X-ray contours (dotted black line) 
are overlayed. The Radio contours start at 4(t (Icr = 50 AtJy). Consecutive Radio and X-ray contours double in value. Surface brightness 
maps for all detected lines are shown in Appendix |X] The northern field is not shown. 
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Figure 7. ABELL 2597 Velocity. We show velocity maps for the H2 1-0 S(3) (top) and Pa a (bottom) lines in units of km s~^ with 
respect to the systemic velocity of the BCG. The velocity of the gas is derived by fitting the spectrum for each spatial pixel in the data 
cube by a single Gaussian. The data was smoothed by two pixels in the spatial plane. The stellar nucleus is indicated by the black cross. 
VLA 8.4 GHz Radio Continuum contours (solid black line. ISarazin et al.lll995l ) and CHANDRA X-ray contours (dotted black line) are 
overlayed. The Radio contours start at Act {la = 50 /^Jy). Consecutive Radio and X-ray contours double in value. Velocity maps were 
made for all detected emission lines. These all show the same structure and are thus not all shown here. The northern field is not shown. 
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Figure 8. ABELL 2597 Velocity dispersion. We show velocity dispersion maps for the H2 1-0 S(3) (top) and the Pa a (bottom) lines 
in units of km s^^. The velocity dispersion ctv of the gas is derived by fitting the spectrum for each spatial pixel in the data cube by 
a single Gaussian. The data cube was smoothed by two pi xels in the spatial p lane. The stellar nucleus is indicated by the black cross. 
VLA 8.4 GHz Radio Continuum contours (solid black line. ISarazin et allllQQSl) and CHANDRA X-ray contours (dotted black line) are 
overlayed on the dispersion maps. The Radio contours start at 4(t (Icr = 50 /^Jy). Consecutive Radio and X-ray contours double in 
value. Gas dispersion maps were made for all detected emission lines. All show the same structure and are thus not all shown here. The 
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Figure 9. SERSIC 159-03 Surface brightness. We show surface brightness maps for the H2 1-0 S(3) (top) and Pa a (bottom) lines in 
units of 10^^^ erg s~^ cm~^ arcsec"'^. The maps are obtained by fitting the spectrum for each spatial pixel in the data cube by a single 
Gaussian. The data cube was smoothed by four pixels in both the spatial and spectral planes. The stellar nucleus is indicated by the 
black cross. VLA 8.4 GHz Radio Gontinuum contours (solid black line) and CHANDRA X-ray contours (dotted black line) are overlayed. 
The Radio contours start at 4fT {la = 25 /tJy). Consecutive Radio and X-ray contours double in value. Surface brightness maps for all 
detected lines are shown in Appendix iBl 
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Figure 10. SERSIC 159-03 Velocity. We sliow velocity maps for the H2 1-0 S(3) (top) and Pa a (bottom) lines in units of km s^^ with 
respect to the systemic velocity of the BCG. The velocity of the gas is derived by fitting the spectrum for each spatial pixel in the data 
cube by a single Gaussian. The data was smoothed by two pixels in the spatial plane. The stellar nucleus is indicated by the black cross. 
VLA 8.4 GHz Radio Continuum contours (solid black line) and CHANDRA X-ray contours (dotted black line) are overlayed. The Radio 
contours start at 4(T (Icr = 25 /^Jy). Consecutive Radio and X-ray contours double in value. Velocity maps were made for all detected 
emission lines. These all show the same structure and are thus not all shown here. 
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Figure 11. SERSIC 159-03 Velocity dispersion. We sliow velocity dispersion maps for the H2 1-0 S(3) (top) and the Pa a (bottom) lines 
in units of km s~^. The velocity dispersion av of the gas is derived by fitting the spectrum for each spatial pixel in the data cube by a 
single Gaussian. The data cube was smoothed by two pixels in the spatial plane. The stellar nucleus is indicated by the black cross. VLA 
8.4 GHz Radio Continuum contours (solid black line) and CHANDRA X-ray contours (dotted black line) are overlayed on the dispersion 
maps. The Radio contours start at Acr (Icr = 25 /^Jy). Consecutive Radio and X-ray contours double in value. Gas dispersion maps were 
made for all detected emission lines. All show the same structure and are thus not all shown here. 
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Figure 12. ABELL 2597 Molecular excitation. Excitation diagrams for the H2 lines in regions Al, A2, A3, A4, A5 and A7. The figures 
show the natural logarithm of the normalised line flux ln(F) versus the upper state temperature Tu. The best-fltting LTE model is given 
by the solid black line line. For qualitative comparison reasons we plot three other H2 excitation models. The dotted red line shows 
a low-density UV fluorescence model by Black and van Dishoeck 1987 (their model 14; n = 3x10^ cm""^, T = 100 K, /yv = 1x10^ 
relative to /tot). The dotted blue line shows a high-density UV fluorescence model by Sternberg & Dalgarno 1989 (their model 2D; 
n = Ix 10® cm" -3, T~1000 K, 7uv = 1x10^ relative to /tot). The green line shows the Ferland et al. 2009 cosmic ray model for the 
Perseus cluster. For region A3 the uncertain H2 1-0 S(2) flux value is shown, but it is not used for the fit. 
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Figure 13. ABELL 2597 Filaments. We show from top to bottom the flux, velocity and velocity dispersion along the filaments in A2597. 
The left images show the northern filament and the right images the southern filament. Each point shown represents a l"xl" area 
along a filament. These areas are marked by the green and red squares in Fig. [3] The black points correspond to green squares and the 
red points to the red squares. The distance of the points along the northern filament are given in kpc from the nucleus, here north is 
positive. The distance of the points along the southern filament are given in kpc from the north-eastern tip of this filament. The black 
points trace the southern filament from its north-eastern tip towards the south-west. The red points trace the southern filament from its 
north-eastern tip towards the south-east. The plus, asterisk, diamond and triangle symbols indicate values obtained for the H2 1-0 S(l), 
1-0 S(3), 1-0 S(5) and the Pa a lines respectively. The blue lines and symbols indicate values obtained for the X-ray emission as observed 
by CHANDRA. These X-ray points have been normalised with respect to the maximum value of the H2 1-0 S(3) line. 
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Figure 14. SERSIC 159-03 Filaments. We show from top to bottom the flux, velocity and velocity dispersion along the filaments in 
S159. The left images show the northern filament and the right images the western filament. Each point shown represents a l"xl" area 
along a filament. These areas are marked by the green squares in Fig. |3] The distance of the points along the northern filament are given 
in kpc to the north-eastern tip of this filament. The distance of the points along the western filament are given in kpc from the nucleus. 
The plus, asterisk and triangle symbols indicate values obtained for the H2 1-0 S(l), 1-0 S(3) and the Pa a lines respectively. The blue 
lines and symbols indicate values obtained for the X-ray emission as observed by CHANDRA. These X-ray points have been normalised 
with respect to the maximum value of the H2 1-0 S(3) line. 
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Figure 15. ABELL 2597 Position- Velocity Diagram. Each point shown represents a 0.25" X 1.0" area along a pseudo long slit with a 
PA of 105.5 degrees. The slit centre is placed 1 kpc south of the nucleus. This is done to avoid the velocity features just north-east 
and north-west of the nucleus. The plus, asterisk, diamond, triangle and square symbols indicate values obtained for the H2 1-0 S(l), 
1-0 S(3), 1-0 S(5), Pa a and Fe II lines respectively. The molecular line velocities are shown in red and the ionised line velocities are 
shown in blue. 
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Figure 16. ABELL 2597 Low velocity shock excitation for region Al. The figure shows the natural logarithm of the normalised line flux 
ln(F), normalisation with respect to the H2 1-0 S(l) line, versus the upper state temperature Tu. The best-fitting LTE model is given by 
the solid black line line. The best-fitting low velocity shock model (shock velocity = 22 km s~^, pre-shock density njj = 5x10'' cm~'^ 
and magnetic field strength B = 5.6 mG) is shown by the purple solid line. The blue dash-dot lines show what happens to the best-fitting 
model if we change the pre-shock density. Increasing the pre-shock density from 1 X 10^ cm' -3, IxlO'* cm-3 to 1x10^ cm ^ the model 
becomes progressively flatter. Similarly increasing the shock velocity from 15 km s~^, 20 km s~^ to 30 km s~^, green dashed lines, flattens 
the model. However, increasing the magnetic field strength from 1.1 mG, 3.4 mG to 7.8 mG, red dotted lines, leads to a progressive 
steepening of the model. 
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Figure Al. ABELL 2597 Surface Brightness Maps. We show surface brightness maps for all detected lines that could be mapped on 
a pixel to pixel basis. The surface brightness is given in units of 10~^^ erg cm~^ arcsec"^ and contours are drawn starting at 
1.5x10"^'^ erg s~^ cm~^ arcsec"^ in steps of 2^/^ with n=0,l,2,... The maps are obtained by fitting the spectrum for each spatial pixel 
in the data cube by a single Gaussian. The data was smoothed by four pixels in both the spatial and spectral planes. Maps are shown 
for H2 1-0 S(0), 1-0 S(l), 1-0 S(2), 1-0 S(3), 1-0 S(4), 1-0 S(5), 2-1 S(3), Pa a, Br 5 and Fe II (1.81 ^lm). The Br 5 map includes the flux 
due to H2 2-1 S(5). The Br 7, H2 2-1 S(2) and H2 2-1 S(4) lines are too weak and are thus not shown here. The northern-most observed 
field is not shown here. 



38 J. B. R. Oonk et al. 




Figure Al. continued 
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Figure Bl. SERSIC 159-03 Surface Brightness Maps. We show surface brightness maps for all detected lines that could be mapped 
on a pixel to pixel basis. The surface brightness is given in units of 10~^^ erg s~^ cm~'^ arcsec"^ and contours are drawn starting at 
1.5x10"^'^ erg s~^ cm~^ arcsec"^ in steps of 2"/-^ with n=0,l,2,... The maps are obtained by fitting the spectrum for each spatial pixel 
in the data cube by a single Gaussian. The data was smoothed by four pixels in both the spatial and spectral planes. Maps are shown 
for H2 1-0 S(l) 1-0 S(3) and Pa a. 
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Figure CI. ABELL 2597 Line Spectra (Region Al). For the Br 7 and H2 2-1 S(2) lines we used a four pixel spatial and spectral 
smoothing. For all other lines a two pixel spatial smoothing was used. The solid black line shows the measured line spectrum. The dotted 
blue line shows the Gaussian fit to the line profile. The dashed red line shows the spectrum from the corresponding spatial median region 
(this spectrum has been scaled to match the area of the selected region). 




Figure C2. ABELL 2597 Line Spectra (Region A2). For H2 1-0 S(0), 2-1 S(3), 2-1 8(2), Br 7, Br <5 wc used a four pixel spatial and 
spectral smoothing. For all other lines a two pixel spatial smoothing was used. The lines and symbols used are the same as in Fig. ICll 
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Figure C3. ABELL 2597 Line Spectra (Region A3). For all lines we used a two pixel spatial smoothing. The lines and symbols used 
are the same as in Fig. ICll 
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Figure C4. ABELL 2597 Line Spectra (Region A4). For H2 1-0 S(2) and 1-0 S(4) we used a four pixel spatial and spectral smoothing. 
For all other lines a two pixel spatial smoothing is used. The lines and symbols used are the same as in Fig. ICll 
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Figure C5. ABELL 2597 Line Spectra (Region A5). For all lines we used a four pixel spatial and spectral smoothing. The lines and 
symbols used are the same as in Fig. ICll 
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Figure C6. ABELL 2597 Line Spectra (Region A6). For all lines we used a four pixel spatial and spectral smoothing. The lines and 
symbols used are the same as in Fig. ICll 
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Figure C7. ABELL 2597 Line Spectra (Region A7). For the H2 1-0 S(3) we used a two pixel spatial smoothing. For all other lines a 
four pixel spatial and spectral smoothing was used. The lines and symbols used are the same as in Fig. ICll 
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Figure C8. SERSIC 159-03 Line Spectra (Region Bl). For all lines a two pixel spatial and spectral smoothing was used. The solid black 
line shows the measured line spectrum. The dotted blue line shows the Gaussian fit to the line profile. The dashed red line shows the 
spectrum from the corresponding spatial median region (this spectrum has been scaled to match the area of the selected region). 
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Figure C9. SERSIC 159-03 Line Spectra (Region B2). For all lines a two pixel spatial and spectral smoothing was used. The lines and 
symbols used are the same as in Fig. IC8I 



H2 1 -0 S(l) 




H2 1 -0 S(3) 



Po alpha 




2,235 2,240 2,245 2,250 



2,060 2,055 2,070 2,075 




,975 1,980 1,985 1,990 



Figure CIO. SERSIC 159-03 Line Spectra (Region B3). For all lines a two pixel spatial and spectral smoothing was used. The lines and 
symbols used are the same as in Fig. IC8I 



44 J. B. R. Oonk et al. 



H2 1 -0 S(1) 




H2 1 -0 5(3) 



Pa alpha 



2.235 2.24D 2.245 2.250 




2.060 2.065 2.070 2.075 



Figure Cll. SERSIC 159-03 Line Spectra (Region B4). For all lines a two pixel spatial and spectral smoothing was used. The lines and 
symbols used are the same as in Fig. IC8I 
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Figure C12. SERSIC 159-03 Line Spectra (Region B5). For all lines a two pixel spatial and spectral smoothing was used. The lines and 
symbols used arc the same as in Fig. IC8I 




Figure C14. SERSIC 159-03 Line Spectra (Region B7). For all lines a two pixel spatial and spectral smoothing was used. The lines and 
symbols used arc the same as in Fig. IC8I 



